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“oleh pA — Ce mémoire a pour but l'étude théorique de la diffusion 
«-18 changement de longueur d’onde de la lumiére par un fluide monoatomique. 
“oix dune méthode statistique:la théorie cinétique, telle qu’elle a été 


codifiée par Gibbs, parait convenir pour traiter les problémes relatifs aux 
liquides. 
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Cuapitre I. — Fluctuations en densité. Ces fluctuations sont lune des causes 
essentielles de la diffusion. — 1. Le calcul des moyennes dans un systéme 
constitué par un grand nombre de molécules ; premier exemple : définition 
de la densité locale, valable aussi bien pour un liquide ou elle est continue 
que pour un cristal ow elle est quasi-discontinue. — 2. La définition des fluc- 
tuations en densité dans un liquide pur ; leur calcul classique ; formule dite 
de Smoluchowski pour le calcul du carré moyen An? dans le cas d’un fluide 
homogéne. — 3. Sur une propriété générale des fluctuations simultanées 
dans deux grands volumes non contigus. — 4. Le probléme de l'état critique : 
la formule de Smoluchowski perd toute signification au point critique. 
Critique des théories proposées jusqu’a ce jour pour résoudre cette difficulté. 
— 5. Des expériences de Gouy donnent la clé du probléme : au point cri- 
tique, le fluide n’est pas homogéne a cause de I’action de la pesanteur ; 
c’est pourquoi la formule de Smoluchowski est défaillante. Le calcul des 
fluctuations dans un fluide rendu inhomogéne par la pesanteur doit étre 
entrepris ; il pose des problémes de mécanique statistique a résoudre d’abord. 
— 6. Les densités simultanées, notion due a Ornstein et Zernike,fondamen- 
tales dans la statistique des milieux moléculaires et particulierement des li- 
quides, sont définies. — 7. Expression du carré moyen de la fluctuation 
en fonction de la densité simultanée double. Expression générale de la densité 
simultanée double en fonction du nombre total des molécules. — 8. Etude 
statistique de la répartition d’un fluide en équilibre thermodynamique dans 
un champ de forces : application de la méthode générale d’Ursell. La for- 
mule (90) donne la densité en fonction du champ ; la formule (94) donne le 
champ en fonction de la densité : elle est trés importante et générale, convient 
pour un gaz, un liquide, pour la zone de séparation liquide-gaz et pour un 
solide cristallisé. L’équation d’état, déja donnée par Ursell, s’en déduit. 
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aux petites distances, (101) valable aux grandes. — 9. Expression du carré 
moyen des fluctuations améliorée d’aprés (101) ; le volume ou sont calculées 
les fluctuations est beaucoup plus petit que le volume total. — 10. Le champ 
de forces est supposé uniforme : c’est le champ de pesanteur ; le carré moyen 
cherché s’évalue alors de maniére compléte a partir de paramétres macros- 
copiques ; le résultat est, satisfaisant méme pour l’état critique. — 11. Nou- 
velle démonstration de l’importante formule d’Ornstein et Zernike, qu 
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exprime le coefficient de compressibilité en fonction de la densité simul- 


tanée. 
Cuapitre II. — Fluctuations en densité et en concentration dans les mélanges 
binaires. — 1. Définition des fluctuations. Leur calcul pour un méiange 


homogéne. Les calculs antérieurs d’Hinstein étaient incomplets, ceux de 
Raman et Ramanathan inexacts. Résultats essentiels formulés en (144), 
(145), (1446), (147), (148). — 2. Au point critique, de miscibilité complete, 
les résultats précédents perdent leur signification. Une théorie, similaire a 
celle développée pour un fluide pur et qui consiste a tenir compte du réle 
de la pesanteur, conduit 4 une évaluation partielle des fluctuations, satis- 
faisante dans tous les cas. Remarques expérimentales. 


Cuapitre III. — Esquisse d’une représentation de l’agitation thermiave 
dans un fluide sous forme ondulatoire qui tient compte, contrairement aux 
conceptions d’Hinstein, de la structure moléculaire des fluides. 
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d’une molécule complexe. —-2. Formule générale (200) de polarisation d’un 
systéme de molécules monoatomiques non polaires dans un champ électrique. 
Théorie générale des milieux polarisés fluides, résumée par (207), (208) 
et (211). Cas d’un fluide homogéne isotrope et d’un champ uniforme : la 
constante diélectrique (220). — 3. Rappel des hypothéses classiques sur la 
polarisation d’une molécule simple placée dans une onde lumineuse et sur 
Ponde diffusée. — 4. Généralisation des résultats de lélectrostatique : 
théorie de la propagation de la lumiére et de la réfraction dans un fluide 
monoatomique transparent, résumée par (236), (237), (238), (241). — 5. Re- 
marques sur la formule de Lorenz-Lorentz. Le développement numérique 
de la théorie montre pourquoi la formule de Lorenz-Lorentz est bien vérifiée 
par les fluides monoatomiques. 
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de la réfraction qui conduit la lumiére vers le milieu d’observation (vide) 
est mis en évidence. — 6. Rapport de Rayleigh. Sa définition. Ses diffé- 
rentes expressions. I] est peut-étre permis d’extrapoler la formule obtenue 
pour des liquides 4 molécules anisotropes. — 7. Résultat essentiel pour le 
rapport de Rayleigh ; résumé des hypothéses ; rappel des notations ; dis- 
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INTRODUCTION 


| bees théorique de la diffusion sans changement de longueur 
d’onde de la lumiére par un liquide en équilibre thermodyna- 
mique n’a été abordée sérieusement jusqu’a présent que par Rocard. 
Cet auteur a montré clairement [1, p. 188 et 2, p. 265] que le pro- 
bléme de la diffusion est étroitement lié a celui de la structure des 
liquides, mais il s’est contenté d’adopter a propos de cette struc- 
ture des hypothéses un peu sommaires si bien que la solution qu’il 
a donnée dans des mémoires d’une lecture attachante n’est pas 
dé finitive. 

Je me propose dans le présent travail de reprendre cette étude. 
Pour faciliter ma tache je me suis limité au point de vue classique 
et j’ai choisi le modéle moléculaire le plus simple possible : les mo- 
lécules sont de forme sphérique et isotropes au point de vue 
optique. Par contre le réle de la structure complexe des liquides a 
été traité de maniére rigoureuse. I] s’est trouvé que le probléme 
de la diffusion, ainsi restreint, aboutit 4 une solution concise et 
adaptée aux calculs numériques. La vérification de la formule 
obtenue relativement a l’intensité de la lumiére diffusée par un 
liquide doit done servir de critérium au sujet des hypothéses que 
j'ai adoptées sur la structure des liquides ; malheureusement les 
liquides dont les molécules réalisent le modéle moléculaire choisi 
et, qui sont les liquides monoatomiques n’ont pas été étudiés expé- 
rimentalement jusqu’a présent : la vérification ne peut donc étre 
tentée que médiocrement 4 l’aide de données relatives a des 
liquides constitués de molécules complexes. Avant d’arriver a ce 
résultat insuffisant il m’a fallu d’ailleurs revoir ou compléter diffé- 
rents raisonnements de mécanique statistique ou d’optique mole- 
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culaire, ce qui a conduit 4 quelques conclusions intermédiaires qui 
sont nouvelles. 

En fait d’hypothéses sur la structure des liquides, j'ai admis 
simplement que les propriétés des liquides en équilibre thermody- 
namique ne sont rien d’autre que des conséquences de la théorie 
cinétique des gaz classique, du moins telle qu’elle a été codifiée 
par Gibbs dans « Elementary principles in statistical Mechanics » 
[3]. Ces principes constituent un corps de doctrine presque tout a 
fait cohérent, auquel on peut reprocher seulement de ne pouvoir 
montrer clairement pourquoi un systeéme matériel quelconque 
abandonné a lui-méme tend vers l’équilibre thermodynamique. 
Cette critique a peu d’importance pour les recherches qui suivent 
puisqu’elles utilisént uniquement des raisonnements relatifs a |’é- 
quilibre. 

Slater [4] a présenté la transposition ondulatoire de la mécanique 
statistique classique ; c’est avec cette statistique seulement que 
pourra s’édifier une théorie tout 4 fait définitive des liquides. La 
statistique classique, beaucoup plus maniable, constitue toutefois 
une bonne approximation lorsque la température n’est pas trop 
basse, parce que les interactions moléculaires ne sont pas trés 
intenses et parce que, avec le modéle moléculaire choisi les rota- 
tions des molécules sur elles-mémes sont exclues. I] ne sera peut- 
étre pas inutile d’ailleurs d’avoir reconnu les phénoménes que peut 
interpréter la statistique classique et dont l’aspect essentiel par 
suite ne sera pas 4 attribuer a quelque particularité encore mal 
connue de la mécanique ondulatoire. Enfin certaines formules 
difficiles 4 démontrer mais élégantes de statistique classique qu’il a 
fallu établir mettront peut-étre sur la voie de formules analogues 
dans la nouvelle statistique. 

La théorie des Jiquides n’est qu’une branche de la théorie géné- 
rale des différents états de la matiére, qui concerne aussi l'état 
gazeux et les états solides cristallin et vitreux. La mécanique sta- 
tistique a servi avec succés a la théorie des gaz : cela a constitué 
la théorie cinétique des gaz ; les chapitres que je vais développer 
de la théorie des liquides méritent donc d’étre considérés comme 
des chapitres de la « théorie cinétique des liquides ». Comme la 
théorie des gaz est beaucoup plus avancée que celle des liquides, 
"habitude s’est répandue d’appeler la théorie des gaz du terme 
trop bref de « théorie cinétique » et aussi s’est répandue l’erreur 
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plus grave de généraliser les résultats de la théorie cinétique des 
gaz et de les appliquer indiment aux liquides. Le présent travail 
a la prétention d’échapper a cette critique. 

Puisque la mécanique statistique réussit avec les liquides, il 
parait vraisemblable qu’elle doit réussir aussi 4 propos de l'état 
solide, en particulier avec l'état cristallisé ; jusqu’a présent pour- 
tant il n’y a pas grand rapport entre la théorie des liquides et celle 
de l'état cristallin ; cela tient & une disposition essentielle de la 
théorie des solides : cette théorie construit d’abord le réseau cris- 
tallin, si bien décrit par les spectres de rayons X, en négligeant 
agitation thermique, ce qui supprime le probléme statistique ini- 
tial posé par les autres états de la matiére ; elle se contente ensuite 
d’envisager les perturbations de petite amplitude qui écartent a 
chaque instant les particules de leur position moyenne. C’est pour- 
quoi la statistique du milieu cristallin a actuellement un aspect 
fonciérement différent de celle des liquides. Les deux points de vue 
devront néanmoins se rapprocher dans la suite : il m’apparait 
comme certain déja que les équations générales que j’établirai sur 
la répartition des molécules dans un liquide s’appliquent sans 
aucune modification au probléme du solide. 

Comme les recherches expérimentales récentes ont montré une 
certaine continuité dans les propriétés d’une substance lors de la 
fusion (voir Magat [5] et aussi Fowler [6]), plusieurs auteurs ont eu 
Pidée non pas de rapprocher la statistique du solide de celle qui 
convient aux liquides mais au contraire d’appliquer aux liquides 
les procédés qui ont réussi avec les solides cristallisés. Comme la 
théorie des solides est plus avancée que celle des liquides, on pou- 
vait espérer résoudre dans cette voie, rapidement et par simple 
analogie, des problemes ardus. (La continuité en question n’est 
peut-étre pas d’ailleurs une propriété générale : c’est du moins 
lopinion de Randall aprés analyse du bismuth aux rayons X[87)}). 

C’est ainsi qu’a pris naissance l’hypothése d’une structure micro- 
cristalline des liquides ; Miiller [7] a par exemple développé récem- 
ment une théorie de l’optique des fluides transparents fondée en 
particulier sur cette hypothése de Stewart [86] : on suppose l’exis- 
tence dans un liquide, autour d’une molécule donnée, d’un groupe 
de 20 4 30 molécules distribuées suivant un réseau ; une molécule 
qui appartient au groupe ne peut occuper qu’une série discontinue 
de positions ; au contraire la position d’une molécule qui n’appar- 
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tient pas au groupe est a peu prés indépendante de la position de la 
molécule centrale. 

Il est difficile de comprendre comment se juxtaposent les diffé- 
rents groupes microcristallins dont l’ensemble constitue le liquide. 
Miiller ne donne aucune précision a ce sujet. Cet auteur a sans doute 
pris Pidée de structure microcristalline trop a la lettre : il y a certes 
une certaine régularité dans la distribution des molécules d’un 
liquide, qui est mise en évidence par I’étude de la diffusion des 
rayons X ; mais cette régularité n’est pas statique comme pour un 
solide, elle n’est que statistique. I] ne semble pas quil y ait avan- 
tage a se représenter cette régularité en déformant l’image que l’on 
a de la régularité cristalline, comme ont essayé de le faire égale- 
ment Prins et Petersen [8] et Bernal [78]. 

L’hypothése microcristalline parait avoir trouvé un renfort dans 
la notion de groupe cybotactique introduite autrefois a la suite de 
étude du point critique : il est certain que dans le domaine cri- 
tique les molécules ont tendance a se grouper les unes au voisinage 
des autres pour former des essaims instables qui laissent entre eux 
un vide relatif. Mais les essaims cybotactiques ne présentent aucun 
caractére microcristallin. D’un liquide ordinaire il est permis de 
dire qwil constitue tout entier un seul groupe cybotactique. 

Une autre conception empruntée a la théorie des solides est celle 
d’ondes d’agitation thermique. Le chapitre III est destiné 4 mon- 
trer comment cette conception doit étre envisagée pour étre incor- 
porée de maniére cohérente dans la théorie moléculaire et déductive 
exposée ici. 

Rappelons enfin les calculs de Bernal et Fowler et de Magat sur 
Peau, de London sur I’hélium au zéro absolu ; les hypothéses d’ An- 
drade a propos de la viscosité, les relations empiriques trouvées par 
Bauer, Magat et Surdin, qui fourniront a la théorie des liquides de 
nouveaux problémes a résoudre, mais qui ne peuvent nous servir 
de point de départ [88]. 

Pour ce qui regarde l’action d’une onde lumineuse sur une molé- 
cule, nous adopterons également des hypothéses trés simples. 
Quoique la théorie ondulatoire de ce phénoméne, qui nous intéresse 
seulement loin de la résonance, soit beaucoup plus avancée que 
celle de Péquilibre thermodynamique — le lecteur se reportera par 
exemple a l’article de Placzek du Handbuch der Radiologie [9] — 
je me suis contenté a ce sujet des idées classiques émises autrefois 
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par Lorentz [10] parce que les deux doctrines ne présentent de 
différences pratiques essentielles que pour le smolécules complexes, 
qui sont exclues a priori des problémes thermodynamiques que 
nous sommes en état d’aborder. 

Il y ala une maniére de procéder qui n’est encore que provisoire : 
en effet le modéle moléculaire sur lequel on a l’habitude de raison- 
ner est extrémement artificiel ; qu’on s’en fasse l'image caricatu- 
rale suivante : une coquille qui posséde des propriétés uniquement 
mécaniques, transparente comme le vide et qui loge en son centre 
un résonateur optique incapable d’actions mécaniques. En réalité 
ce sont les mémes résonateurs périphériques qui sont responsables 
a la fois des phénoménes optiques et des phénoménes mécaniques : 
seule la mécanique ondulatoire est capable de réunir dans la méme 
synthése ces deux classes de phénoménes. 


CHAPITRE I 


Depuis les travaux de Smoluchowski on sait que les fluctuations 
en densité sont l’une des causes de la diffusion de la lumiére dans 
les fluides et la cause essentielle de l’opalescence critique. C’est 
pourquoi je commence par étudier les fluctuations en densité, sans 
faire aucune allusion a l’optique. Rappelons seulement que ce qui 
importe surtout pour l’étude de la diffusion est le calcul du carré 
moyen de ces fluctuations. 


1. Le calcul des moyennes. — Appliquons donc la statistique 
dont Gibbs a donné le modéle ; j’ai déja indiqué la maniére de s’y 
prendre dans un mémoire antérieur [11] ; je renverrai a ce mémoire 
le moins possible parce que je désire développer ici une nouvelle 
méthode d’approximations successives, méthode inaugurée par 
Ursell et qui nécessite un exposé complet. 

Le fluide théorique est constitué par des molécules trés nom- 
breuses, supposées pour le moment toutes identiques, et numérotées 
de1a WV. AVinstanttla molécule numérotée J est caractérisée par les 
coordonnées d’espace et de vitesse 2, y; z; et u; ¥; w, de son centre 
M;. Nous poserons pour abréger 


dw, = dx;dy;dz; dy; = du;dy,dw; 
et aussi 
(1) dws, => dwsdu, dys«c = dy;dyx. 
Un systeme de valeurs de x, y, z et de u, », w donné pour un 
certain instant ¢ définit une phase du fluide. A chaque phase pos- 
sible correspond un coefficient de probabilité D,fonction des 2N + 1 


variables précédentes ; D est une fonction continue, dérivable, 


symétrique par rapport aux diverses molécules, satisfaisant a 
Pé quation 


(2) J [ Daerendason == Be 
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de plus une grandeur F quelconque qui est fonction de la phase, 
autrement dit des 2N + 1 variables, a une valeur moyenne F 
fonction du temps donnée par la formule 


(3) F = J [ DF avraadzstn 3 


seules les valeurs moyennes sont accessibles a l’expérience. II est 
entendu que dans les intégrales précédentes les vitesses varient de 
—o a-+ o et que le centre de chaque molécule prend toutes les 
positions 4 l’intérieur du récipient qui contient le fluide. 

Comme application de ces notions définissons la densité moyenne 
dans une fraction A quelconque du volume total V du récipient ; 
considérons la fonction 


F(x, yz, Z3) ou plus simplement F(M;) ou Fy; 


qui est égale 4 1 lorsque le point M, se trouve dans A et qui est 
égale a 0 lorsqu’ilest al’extérieur ; le nombre de molécules qui se 
trouve dans A est 


(4) > F. 
1 


Calculons sa valeur moyenne n ; on a d’abord évidemment 


(5) F, = Hi ff Ddwydtng3...5712--n, 


les limites d’intégration sont les mémes que tout a I’heure saul 
pour l’intégration qui concerne l’élément dw, et qui est limitée au 
volume A ; par raison de symétrie, on a 


et 
(6) n = NF}. 
Posons alors 


(7) y= N ff Ddowndia2--n' 


les limites d’intégrations sont toujours les mémes que pour les 
intégrales (2) ou (3) (*) ; 4 ne dépend que des coordonnées 


71, Y1; 21; 





(*) Lorsqu’il en sera ainsi par la suite, les limites d’intégrations ne seront 
pas explicitées en général. 
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et peut s’écrire d’une maniére plus explicite 

va = V(L1, Yr, 21) 5 
cette fonction recoit le nom de « densité moléculaire au point M, ». 
Il est clair que l’on a d’aprés (2) 


(8) [onde = N 


et d’aprés (6) 


(9) n= [ vide. 


vA 


2. La définition des fluctuations et leur calcul classique. — 
Avant de les perfectionner, il faut reprendre les calculs des fluc- 
tuations sous leur forme classique, ce qui va nous obliger 4 évaluer 
les moyennes d’une maniére moins simple que par (3). Nous ne 
reviendrons 4 lemploi de (3) qu’au paragraphe (6). 

Le volume A contient done un nombre de molécules 


=F; 
variable d’une phase a l’autre ; posons 
(10) An = =F;—n; 


An est par définition la fluctuation du nombre des molécules con- 
tenues dans A pour la phase considérée. II est clair que cette fluc- 
tuation est nulle si le volume A se confond avec le volume V et que 
de toute maniére 


(11) An = 0. 


La fluctuation du nombre des molécules est accompagnée d’une 
fluctuation de la densité moléculaire 


An/A; 


ces deux fluctuations sont dans un rapport constant ; étudier l’une 
c’est étudier l’autre. I] est dans Pusage de considérer la seconde de 
ces grandeurs plutét que la premiére. J’ai trouvé plus simple au 
contraire de raisonner en général sur la fluctuation du nombre des 
molécules qui sera souvent désignée ici, pour abréger, du seul 
terme de fluctuation. 

Recherchons quelle est la probabilité d’avoir a un instant donné 
n + An molécules dans A et par suite VN —n—An dans le volume 
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restant V — A ; pour la symétrie des notations posons temporaire- 
ment 


M=n-+An ng = N—n—An B=V—A; 


imaginons d’abord que les n, molécules en question se trouvent 
étre justement les molécules numérotées depuis 1 jusqu’a n, ; la 
probabilité d’avoir cette répartition est donnée par l’intégrale 


(12) 2} =|. . fof + Ae -[, [ Daordos . ‘dion, -din, > . dwyd¥ 19.008 3 
A A B B 


supposons maintenant que les n, molécules contenues dans A ont 
des numéros quelconques ; il y a 


N! 
nal Ng ! 
maniéres de répartir V molécules entre deux récipients A et B de 
fagon que n, quelconques d’entre elles soient dans le récipient A 
et le reste dans B. La probabilité cherchée est donc 
Ni 


Nya! ng! 


(13) 


Mais la probabilité d’une fluctuation ne nous intéresse jusqu’a pré- 
sent que dans le cas d’un fluide en équilibre thermodynamique ou 
pour préciser d’un fluide en équilibre dans un thermostat (voir a ce 
sujet L. Brillouin [12]) ; on sait alors depuis Gibbs que le coefficient 
de probabilité D est donné par |’expression 


(14) D= el¥ — 2)/kT 


ou £ est l’énergie totale du systéme de molécules considérée comme 
une fonction des coordonnées de vitesse et de position des molé- 
cules,7 la température thermodynamique absolue, k la constante 
de Boltzmann * et ¥ le potentiel thermodynamique 4 volume 
constant qui figure comme une constante dans les calculs de 
moyennes et qui est défini par la condition (2). 

Supposons maintenant que le volume A, soit comme le volume 
V lui-méme, grand devant les dimensions moléculaires ; supposons 
en outre que la surface qui sépare les volumes A et B soit a peu 
prés lisse ; alors l’énergie EL se compose de trois termes: £, énergie 


———— NI 
(*) k est encore le rapport R/No, R constante des gaz parfaits,Vonombre 
d’Avogadro. 
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des molécules contenues dans A, E, énergie des molécules con- 
tenues dans B, Ey énergie mutuelle des deux groupes de molécules, 
dont le dernier est négligeable devant les deux autres. On peut 


donc poser 
(15) E= Ey -+- Es, 
l’intégrale (12) se sépare en un produit de deux intégrales et 
donne 
is — CHT fo FA Aopen dato 
(16) 


] ff mF dem roading ; 
B 


mais le volume A considéré isolément avec ses n, molécules a un 
potentiel thermodynamique 4, défini par une condition équiva- 
lente a (2) qui s’écrit 


(17) ff ela Balk Pdrendzizom aE 


ou 
(18) e Palkt = [femdom drat 


il y ade méme pour B un potentiel thermodynamique ¢, ; Ys est 
fonction de n, et de T comme ¢, est fonction de n, et de T. La 
probabilité (13) s’écrit donc 


(19) Ni el¥—Ya—da/kT yy Nal: 
a laide de la formule de Stirling on remplace 


nlogn—n, 
n! par e? °8 . 


on pose 
(20). Y=¢—nlogntn 
et la probabilité s’écrit, 4 un facteur constant prés, 
(24) exp — (Wy + Wp)/kT ; 


¥ comme v est un potentiel thermodynamique ; pour les distin- 
guer Gibbs [3, p. 176] qualifiait ) d’exposant relatif aux phases 
« spécifiques » et Y d’exposant relatif aux phases « génériques » ; il 
suffit de dire ici que YY se confond avec le potentiel thermodyna- 
mique & volume constant souvent formulé. 


(22) Tis P79 
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(S entropie) et que dans le cas ou le fluide est homogéne au sens 
thermodynamique ¥° est proportionnel au volume du fluide envi- 
sagé ; VY, peut donc se mettre sous la forme 
(23) VY, = A®(T,n,/A), 
® étant une nouvelle fonction. Une telle propriété n’appartient pas 
a 

Désignons respectivement par n,, ng les valeurs de n, et de np 
qui rendent la probabilité (21) maximum ; ces valeurs vérifient les 
relations 

n, + ng = N 

(24) OW, /on, = a g/ans ; 


au cours des fluctuations les écarts 
ny — n, et Ng — Ng = — (m — n,) 


restent trés petits ou du moins de grandes valeurs de ces écarts 
sont trés peu probables ; il est légitime de remplacer dans (21) ¥, 
et , par des expressions approchées ; par exemple on peut poser 


(25) Wy = Wy + (0W',/amg) (ma— M1) + 5 (2a /2%Q) (Ma— Mg)? Ho 

Nous allons supposer pour simplifier que le volume B est beau- 
coup plus grand que le volume A : B est vis-a-vis de A un immense 
réservoir 4 molécules, la fluctuation dans B est relativement beau- 
coup plus petite que dans A et on peut se contenter de remplacer 
Ys, par 

Ws \-F (a'¥g/ans) (m3 — Nz) 3 

- limitons d’ailleurs (25) aux termes du second degré par rapport a 
la fluctuation ; dans ces conditions la probabilité (13) ou (21) s’é- 
crit 
(26) Ce — (2°, /02n,) (ng — n,)?/2kT ; 
C est un facteur constant. Posons pour abréger 
(27) a = (2W/a*n,)/2kT, 
C est défini par la condition 


(28) ECo— Aa — 7, )* = A, 


dans la somme sous signe © n, prend toutes les valeurs de 1a JN ; 
la valeur moyenne de ny, est alors 


(29) Na =i = ECn,e7 W421" « 
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le carré moyen de la fluctuation est 
(30) An? = (n, — n)? = 2C(n, — n)®e— *(*a— 7)” 

On sait, a priori, et nous nous sommes déja servi de ce fait, que 
les fluctuations un peu importantes sont peu probables, il faut 
done admettre que a est positif et qu'il n’est pas trés petit ; on 
peut remplacer les sommes qui figurent dans (29), (28), (30) par 
des intégrales approchées. (28) devient 


+o , 
(31) ‘ Con aR idap esl 
ou 

4 
(32) C = (r/a)?; 
(29) devient 
+o : 

(33) n=C i ne A —2,) dn. 


et montre que n et nz doivent étre confondus avec les approxi- 
mations que nous avons faites. Enfin (30) devient 


Bik. r++ 
(34) ve ee J (4 — n)2e—2("a—2)*gy 
et donne 
(35) An? = 1/20. @ 


Ce genre de résultat figure déja dans la mécanique statistique 
de Gibbs [3]; je n’ai fait que suivre en résumant le calcul de 
Gibbs ; Smoluchowski [13] a présenté depuis dans son mémoire 
classique les choses d’une maniére un peu différente : au lieu de 
considérer les variations du nombre des molécules dans un volume 
fixe, il a évalué les variations d’un volume contenant un nombre 
invariable de molécules ; le point de vue que j'ai adopté ici est 
Sans doute préférable parce qu'il s’adapte sans modification au 
cas d’un mélange de plusieurs corps purs et si j’ai repris les rai- 
sonnements avec quelques détails, c’est justement afin de pouvoir 
réviser plus loin les affirmations qui ont été publiées a propos des 
mélanges et qui ne me paraissent pas exactes. 

pe reste a expliciter le paramétre a dans un cas important, celui 
ei "3 oe A contient un fluide homogéne ; la pression en 


(36) ow ,/2A = — Pp: 
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en s’aidant de (23) et de la définition générale du coefficient de 
compressibilité 


. 
(37) "yates = (*) 
on trouve que 
(38) a) = 2kTbn2/A 
et que 
(39) An® = kT8n2/A. 


Désignons par p la densité définie comme le nombre de molé- 
cules gramme contenues dans l’unité de volume, nous arrivons a 


(40) An? = RT fen, 


ce qui correspond au résultat de Smoluchowski (celui-ci calculait 
en réalité non pas le carré moyen de la fluctuation,mais la moyenne 
de sa valeur absolue). 


3. Association de volumes. — La formule (40) présente une im- 
portante propriété qui, physiquement, était 4 prévoir sans rai- 
sonnements complexes : 

Imaginons de découper le volume total V en trois parties, 
A, A’, B; A’ est comparable au volume A déja envisagé, B est 
beaucoup plus grand ; nous supposons que A et A’ sont éloignés 
Pun de l’autre et nous n’avons pas 4 admettre pour le mofnent 
que le fluide est thermodynamiquement homogéne. Désignons 
par An et An’ les fluctuations qui se produisent dans A et dans A’, 
respectivement, au méme instant. On a 





(41) (An + An’)? = An? + An 4+ 2Andn’ ; 


mais, puisque A et A’ sont éloignés et « noyés » en quelque sorte 
dans le grand réservoir 4 molécules B, il parait évident a@ priori 
que l’on doit avoir 


(42) AnAn' = 0; 
alors (41) devient 
(43) (An + An’)? = Sn? + On?. 


Pour le cas particulier ou le fluide est thermodynamiquement 





(*) Ici @ n’a bien entendu plus la méme signification que ci-dessus. 
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homogéne, la formule (40) donne aussi bien les fluctuations dans 
A, dans A’, que dans l’ensemble A + A’;les résultats auxquels 
elle conduit pour ces trois volumes vérifient, comme il convient, 
la relation (43). 


4. Le probléme de l'état critique. — Le résultat [40] n’a aucun 
sens pour un fluide a I’état critique puisque pour cet état le coeffi- 
cient de compressibilité 6 est infini. Smoluchowski avait déja 
vu cette difficulté et avait cru la résoudre en utilisant dans le dé- 
veloppement en série (25) jusqu’aux termes du quatriéme degré 
par rapport a ns — na; pour l'état critique on a 
(44) V/an? = 0 ew/ane = 0 
et 

¥, = Vy + (2¥,/on,)(n, —n, ) + (at¥/an,4)(n, — n,)4/24. 

On trouve alors que 

ae: Tr 1/2 
(45) An? = al sop n3/2 ; 


7 est une constante qui a Ja méme valeur pour tous les fluides. 
Ce résultat (45) ne parait pas satisfaisant non plus, cette fois-ci 
parce qu il n’est pas conforme A (43). 

Ce sont Ornstein et Zernike [14] qui se sont rendu compte les 
premiers de ce défaut ; soucieux avant tout d’établir une théorie 
de la diffusion de la lumiére, ils n’ont d’ailleurs pas cherché a 
corriger d’une nouvelle maniére la formule (40), mais seulement 
a lui faire jouer un rdle pas trop génant dans la théorie de la diffu- 
sion. Le probléme précis, qui ne concerne pas optique, de trouver 
une formule de fluctuation dans les grands volumes pour I’état 
critique reste donc entier apres les mémoires d’Ornstein et Zer- 
nike ; il convient d’ajouter que ceux-ci ont eu le mérite considé- 
rable d’attaquer le probléme des fluctuations dans les petits vo- 
fumes avec une méthode, a peu prés, moléculaire et qui a trouvé 
sa plus belle application dans le mémoire fondamental de Zer- 
nike et Prins [15] sur la diffusion des rayons X par les liquides. 

Plus récemment, Rocard [16] a proposé une nouvelle conception 
de état critique, qui se rattache a une théorie générale de l'état 
liquide dont lorigine remonte & un mémoire d’Einstein [17] et 
qui a été développée apres Einstein par L. Brillouin [18]. Voir 
aussi L. Brillouin [19] et Lucas [20}. 
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Cette théorie consiste 4 représenter l’agitation thermique dans 
les liquides par des trains d’ondes acoustiques dont Yamplitude 
moyenne est nulle, mais qui ont en moyenne une certaine énergie. 
Il est certain que cette conception représente a priori une bonne 
part de vérité puisque tout mouvement peut s’analyser sous forme 
d’ondes sinusoidales a l’aide de la série de Fourier. Pour aller 
plus loin, Einstein a supposé que les ondes sont harmoniques, A 
chacune d’elles correspond une énergie cinétique et une énergie 
potentielle égales ; l’énergie totale est la somme des énergies de 
chaque onde. On peut imaginer des ondes libres, pour lesquelles 
énergie potentielle est nulle. La fréquence de ces ondes est limi- 
tée supérieurement de maniére que le nombre de degrés de liberté 
relatif aux ondes soit égal au nombre de degrés de liberté déduit 
de la structure moléculaire. Il est alors facile de procéder a la 
statistique des ondes et de leur appliquer le principe d’équipar- 
tition de l’énergie entre les différents degrés de liberté. Pour ob- 
tenir des résultats valables aux basses températures, il suffit de 
remplacer le principe d’équipartition par la formule de Planck. 

Cette théorie de l'état liquide se rattache étroitement a la théorie 
analogue développée d’abord par Debye pour les solides. Elle 
se heurte A une difficulté essentielle que ne rencontre pas aussi 
nettement la théorie du solide : les longueurs d’onde acoustiques 
que l’on est amené a envisager sont comparables a la distance 
de deux molécules voisines. I] n’y a aucune vraisemblance, dans 
le désordre de |’état liquide, pour que les ondes correspondantes 
soient harmoniques, pour que_leur énergie totale soit la somme de 
leurs énergies partielles. L’application de l’hypothése des ondes 
d’agitation thermique a la théorie de la diffusion de la lumiére 
montre d’ailleurs clairement les limites de validité de cette hypo- 
thése : elle permet de prévoir de maniére satisfaisante (L. Brillouin 
[21]) la structure fine de la lumiére visible diffusée par un liquide, 
mais elle échoue a propos de la diffusion des rayons X. 

Il faut donc considérer les ondes d’agitation thermique d’aprés 
Einstein comme une approximation mathématique commode, 
valable seulement quand la longueur d’onde acoustique n’est pas 
trop petite, qui ne prétend pas a la rigueur et qui doit céder le 
pas aux principes stirs de la mécanique statistique dés que ceux-Cl 
ne conduisent pas 4 des calculs inabordables. C’est pourquoi je ne 


puis admettre le point de vue de Rocard — il faut cependant 


Yvon I. : 
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retenir de son travail la notion, indiscutable théoriquement, de 
forces capillaires. 

D’une maniére plus générale, la théorie des ondes d’agitation 
thermique mérite le reproche de n’étre pas moléculaire ; elle parait 
inapte a interpréter des phénoménes comme la biréfringence élec- 
trique ou magnétique qui dépendent directement de la structure 
moléculaire. Cet inconvénient a été sensible méme dans la théorie 
des solides puisque Born et Karman [22] ont jugé utile de perfec- 
tionner la théorie de la chaleur spécifique des solides de maniére 
A tenir compte de la structure moléculaire de ces corps. Le cha- 
pitre 11 est destiné 4 montrer qu’il est possible de donner des 
ondes d’agitation thermique une définition conforme 4 la struc- 
ture moléculaire des liquides et 4 la mécanique statistique ; ces 
ondes ne présentent pas la méme simplicité intuitive que dans la 
théorie ancienne : j’ai remis a plus tard d’en tirer des applications 
parce que ces ondes ne conduisent 4 aucune conclusion qui ne soit 
déja impliquée dans la statistique générale. 

Quoiqu’il en soit, il nous reste 4 trouver une inexactitude dans 
le raisonnement de Smoluchowski. 


5. Le réle de la pesanteur. — Je crois que cette inexactitude a 
pour cause l’oubli du réle de la pesanteur. 

Reportons-nous d’abord a de belles expériences de Gouy qui 
datent de 1893 [23]. Cet auteur a chauffé avec une extréme len- 
teur — il estimait lélévation de la température 4 0°0001 par 
heure, un grand nombre de tubes de Natterer diversement 
remplis d’anhydride carbonique. 

La théorie élémentaire du phénoméne affirme que les obser- 
vations doivent étre les suivantes : le ménisque qui sépare a basse 
temperature le fluide en deux couches se déplace quand on 
chauffe le tube ; pour les tubes peu remplis de fluide, le ménisque 
descend et finit par disparaitre a la partie inférieure ; pour des 
tubes plus fortement remplis au contraire, le ménisque monte 
et vient disparaitre 4 la partie supérieure ; enfin, pour un rem- 
plissage unique, intermédiaire entre les précédents, le ménisque 
reste stationnaire au milieu du tube et s’y estompe brusquement 
a une certaine température, a ce moment-la, il ne reste plus qu’une 
phase dans le tube : c’est l’état critique. D’aprés ces conceptions, 
il est invraisemblable qu’un expérimentateur puisse observer 
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état critique puisque cette observation nécessiterait : 1° un rem- 
plissage mathématiquement exact ; 2° l’absence de la plus minime 
inégalité de température, car celle-ci aurait pour effet de 
provoquer d’intenses phénoménes de convection qui trouble- 
raient |’équilibre. 

Cette description théorique choque dés l’abord le sens physique 
de la continuité, et en effet Gouy a vu autre chose : la tempéra- 
ture était assez bien définie pour éliminer les convections acci- 
dentelles ; le remplissage 4 réaliser pour observer la disparition 
du ménisque ailleurs qu’aux extrémités du tube n’est pas diffi- 
cile a réussir puisque, avec des tubes de 27 centimétres de hauteur, 
il n’est défini qu’a 7 % prés. La disparition du ménisque n’est 
pas subite 4 une température déterminée : lorsque la température 
s’éléve, le ménisque se change progressivement en une zone qui 
présente le phénoméne de mirage et dont la hauteur peut atteindre 
1 a 2 millimétres ; cette zone s’estompe ensuite peu a peu. Cette 
modification et cette disparition se produisent au total dans un 
intervalle de température a 0°001. 

I] est facile de donner une interprétation de ces phénoménes, 
et Gouy n’y a pas manqué. Cette interprétation repose sur le fait 
omis par la théorie élémentaire que dans un fluide soumis 4 la 
pesanteur, la densité varie avec l’altitude ; ce phénoméne est 
d’ailleurs d’autant plus accentué que le fluide est plus compres- 
sible, ce qui est justement le cas d’un fluide dans le domaine cri- 
tique. Raisonnons sur un tube convenablement rempli et que l’on 
refroidit tres lentement — il est plus commode de _ rai- 
sonner sur le refroidissement. A 1° au-dessus de la température 
critique, le fluide est 4 peu prés homogéne ; 4 0°01 la diminution 
de la densité de bas en haut est appréciable ; ce n’est d’ailleurs 
pas une fonction linéaire de l’altitude : la compressibilité du fluide, 
fonction de sa densité, n’est pas la méme @ tous les niveaux ; a 
un certain niveau C, qui dépend peu de la température, la com- 
pressibilité est maximum : au voisinage de ce niveau le gradient 
de densité est plus grand qu’ailleurs. A 0°001 de la température 
critique, le phénoméne s’accentue, la variation de densité est 
assez rapide au voisinage du niveau C pour produire le mirage. 
Enfin, a la température critique — dont la définition comporte 
une légére imprécision — la zone de mirage est remplacée 
par une couche capillaire qui sépare le fluide en deux phases : 
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en méme temps pour chacune des phases la compressibilité a 
diminué, elles sont A peu prés homogénes. Quant a la condition 
de remplissage, il est facile de voir qu’elle est trés vague ; avec 
un tube de hauteur infinie en particulier, elle se réduit 4 une iné- 
galité. 

Le phénoméne peut se mettre en équations. Désignons par p 
la densité au niveau d’altitude z, par P la pression correspondante, 
par M la masse de la molécule gramme, par g la valeur absolue de 
Yaccélération de la pesanteur. La température est légérement 
supérieure a la température critique. L’hydrostatique donne la 


relation 

(46) dP/dz = — Mgo 
ou 

(47) aP/ae-de/dz = — Mgo 
ou 

(48) de/dz = — Mgo?S, 


Cette derniére équation montre que p¢ varie rapidement avec 
le niveau puisque 6 est grand au voisinage de l’état critique. 
Recherchons le niveau C au voisinage duquel la variation est par- 
ticuliérement rapide ; pour ce niveau on a 


d?o/dz? = 0 

ou d’aprés (47) 

(49) oP/ao — ¢-a*P/ac? = 0; 

sur lisotherme il existe une densité p- pour laquelle on a 
a Plant ==. 0; 


d’autre part, pour cette densité, suivant la forme classique des 
isothermes, OP /dp est trés petit, si bien que cette densité vérifie 
pratiquement (49) : le niveau C est le niveau de la densité pe. 

Ces remarques cessent d’étre valables quand la température 
se rapproche beaucoup de la température critique, c’est- 
a-dire 4 partir du moment ot expérimentalement apparait la 


zone de mirage : dans la zone de mirage, le rapport 
aP/d 


ne peut plus étre déduit du réseau des isothermes ; la pression 
au niveau z ne dépend pas seulement de la densité ¢ 3 i 
p ensite ¢ a ce niveau, 


of FLUCTUATIONS EN DENSITE I-21 


mais aussi des actions capillaires exercées par les couches voisines. 
Pour ce cas, Rocard a indiqué [24], a la suite d’une théorie de la 
capillarité, des équations correctes dont il a remis a plus tard la 
discussion. Je n’ai pas l’intention de faire davantage sur cette 
question parce que les phénoménes paraissent suffisamment clairs, 
mais il reste & ce sujet un probléme mathématique curieux a 
traiter. 

De toutes ces explications, il résulte qu’il n’existe pas de fluide 
qui corresponde a la définition classique de l’état critique, état 
fluide homogéne a coefficient de compressibilité infinie ; par con- 
séquent la formule (38) qui s’applique seulement a un fluide homo- 
géne, et la formule de fluctuation (40) qui en découle, ne doivent 
pas étre appliquées a cet état. Du méme coup, elles ne comportent 
plus aucune objection sérieuse. 

I] reste maintenant 4 trouver une expression du carré moyen 
des fluctuations qui s’applique au voisinage immédiat de |’état 
critique ; c’est 1a un exercice destiné a assouplir les méthodes de 
la mécanique statistique plutét qu’une recherche susceptible de 
vérifications expérimentales : ces fluctuations en densité se tra- 
duisent bien par la diffusion de la lumiére, mais |’on ne passe pas 
assez simplement de l’une a l’autre pour qu’une théorie des 
fluctuations puisse étre considérée comme une théorie de la 
diffusion. 

On peut se demander ce qu’il adviendrait si le champ de pesan- 
teur disparaissait. Il semble 4 premiére vue qu’alors la difficulté 
représentée par (40) reparait. Je ne le pense pourtant pas, car le 
fluide est toujours soumis 4 un champ, trés localisé mais trés 
intense, le champ da Aa la paroi du récipient, qui serait respon- 
sable a l’intérieur du fluide d’inhomogénéités analogues a celles 
dues habituellement 4 la pesanteur. Ce point est repris plus 
loin avec quelques détails. 


6. Les densités simultanées. — Pour calculer les fluctuations en 
densité d’une maniére générale, j’abandonne la méthode du 
§ 2; ce qu’on va lire est la suite du § 1 ; aucune hypothése parti- 
culiére n’est faite pour le moment sur le coefficient de probabi- 
lité D, sauf que je suppose, pour simplifier, qu’il ne dépend pas 
du temps. 

Envisageons deux fractions A et B du volume total V ; elles 
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n’ont aucune partie commune et, 4 part cela, sont quelconques. 
D’autre part, considérons les molécules deux 4 deux : le nombre 
de couples de molécules est égal a 


(50) SMW — 1). 


Proposons-nous de compter le nombre de ces couples qui ont 
une molécule dans A et l’autre dans B: a cet effet, introduisons 
une fonction F, de M, et de M, qui est égale 41 quand M, est 
dans A et quand M, est dans B et qui est nulle dans les autres 
cas ; il est clair que dans tous les cas, le produit 


Fy. y Fa 
est nul. Le nombre de couples cherché est 
(51) ZE(P yo + Foy) ; 


la somme double comporte autant de termes Fy, + Fy. qwilya 
de couples de molécules. La valeur moyenne du nombre de couples 
de molécules est donc 


1 
5 V(N —- ) ff as — Fy) Ddo4p...~d 712... 


ou, puisque toutes les molécules jouent le méme réle, 


(52) N(N = 1) | f FaDaers vd. 


on supprime le facteur F,, cependant que l’extension de dw, est 
limitée a A et celle de dw,» & B; le nombre cherché devient 


(53) N(N —1) [ if if { Dd dergdengs...xd 79.004. 


L’intégrale 
(54) NN —1) J | Ddusgesd:ayen 


ne dépend que des coordonnées Yi%, TYo22 ; désignons-la par 


ULYsZ, XeYo%q) OU v9; le nombre moyen de couples cherché est 
alors 


(55) | f vrsdondos 
VYASB 


J’ai donné & v4». le nom de densité simultanée aux points M, 
et M,. C’est une généralisation de la densité ordinaire »(M) ; 
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elle peut se généraliser 4 son tour et il y a lieu d’envisager des 
densités relatives a 3, 4,... points a la fois qui se notent v4.95, Yio4, +++ 
Par exemple 


(56) Y123 = N(N — 1)(N <te 2) [ [ Dadvemsdyaa 


En rapprochant (54) et (7) on voit que 


(57) ye =| vi2dwg, 
de méme 
(58) (N — 2)v2 = if ¥493dw3, 


(N — 3)y128 = f rusdos 


Dans un fluide homogéne et isotrope 4. ne dépend que de la 
distance 
32 >= M\M2 


que nous désignerons, quand aucune confusion ne sera possible, 
simplement par r. Hormis le cas des gaz peu comprimés, de telles 
densités ont été peu considérées jusqu’a présent. L’idée en revient 
a Ornstein et Zernike [14] : 4 propos des fluctuations, ils ont uti- 
lisé la densité simultanée 4 deux molécules dans le cas d’un fluide 
homogene indéfini ; et Zernike et Prins [15] l’ont appliquée avec 
bonheur a la diffusion des rayons X par les liquides. Je retrou- 
verai, en les généralisant, une partie des résultats d’Ornstein 
et Zernike, mais je serai obligé de critiquer les autres. 


7. Fluctuations en densité et densités simultanées. — Reprenons 
la fonction F, qui a servi a définir la densité simple. Le nombre 
de molécules contenues dans le volume A est, pour une phase 
donnée, la somme 

LF; 


sa valeur moyenne est 


(9) , n = {nds 
A 


la fluctuation est 


(59) An = =F; — [doy 
A 
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La valeur moyenne de la fluctuation est nulle, mais il n’en est 
pas de méme de son carré moyen ; le carré de la fluctuation est 


2 
(60) An? = =F? + 2==F,F, — 2=F;- fondo +[ [ 4 . 
A x 
on remarque que 
(61) F,? = F,, 


2 . 
(62) [ f nde | _ | ‘fi ¥y¥edw,dw, = af ¥4%0049, 
a id arat A 


et que la valeur moyenne de 2=F,F,, somme qui comprend 


5 MN — 1) termes simulaires, peut s’écrire 
(63) N(N — 1) ff Fy FyDde,.. dey...» 


= a FF. 29d). = if ¥ 9h 19; 
4 


le carré moyen de la fluctuation est done, en tenant compte de 
la valeur moyenne (9) de F,, 


(64) An? = { »do, oo [Ou — ¥y%q) dw 49 ; 
A e 


mais d’aprés (8) mis sous la forme 


fade == 1V 
“ 
et d’aprés (57) on a 


(65) f fb _— 449) dw, dw, = — f vide 
AYJV A 


et d’autre part 


i a (%32—%y¥o) dw, dw, = fp (%32—Vy%q) wy. as if af (Y¥y2— ¥y¥q) dn dwg, 
AJ V A AJ V—A4 


donc 


(66) An? — —f if (vig — ¥4¥9) dw, dus ; 
Av (V—A) 


ce dernier résultat montre, comme il convenait, que le carré 
moyen de la fluctuation est le méme pour A et pour V — 4, 
On pourrait croire que, parce que V’influence d’une molécule 
sur une autre diminue trés vite avec l’éloignement, la différence 
“1g — %4%_ tend rapidement vers 0 quand la distance augmente. 
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Ce n’est pas tout 4 fait exact. Nous .dmettons pour le moment 


et nous montrerons plus loin que »,. a pour expression une série 
entiere en 1/N 


(67) Yyq = YyVe(4yq + 5y2/N + Cy9/N? + ---), 

OU Ayo, Dye, Cro, --- Ne dépendent que des points M, et M, et des 
particularités du fluide dans le voisinage de ces points. Pour la 
plupart des problémes, il est permis de négliger tous les termes 
de la série sauf le premier ; ici nous devons conserver aussi le 
second. Quand la distance augmente, a,, et b,, tendent rapidement 
vers une limite ; la limite de a,, est 1, celle de »,. — yr, est done 
(68) v%_D,9/NV. 


(57) s’écrit avec une approximation suffisante en profitant de (8) 


(69) (N—1)y, = wf vale — 1)dw, + Nv, + ws { vebrdogl 


ou 
(70) Hi ede New | ie secetidere: 


le premier membre est fini et comparable au second parce que 
b,, tend vers une valeur non nulle quand la distance augmente. 

Appliquons nos différentes formules au cas d’un fluide homo- 
gene isotrope, négligeons l’effet perturbateur de la parol qui dé- 
truit toujours ’homogénéité dans son voisinage ; v est constant, 
v2 ne dépend que de r; dans |’intégrale du premier membre de 
(70), b:, peut étre regardé comme constant, cela ne constitue une 
erreur qu’aux trés petites distances ; on a alors 


f vabusdos = bie { oad = Nbyp. 
et 


(71) by. =—1 — [ver — 1)derg ; 
(66) devient d’aprés (67) 


N—n 
(72) An? = — af { (ay. — 1)dw, dw, — i bie 
ALVA 


en remarquant que 
1p vey = N— ny; 
V—A 
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eu égard aux changements de la région A, le second terme dépend 
simplement du nombre moyen de molécules qu'elle conUenys 
le premier terme comporte en facteur a,, — 1 qui s’annule trés 
vite, quand la distance augmente, les domaines d’intégration A 
et V — A peuvent étre réduits 4 deux zones trés minces longeant 
la surface limite de A de part et d’autre ; si cette surface n’est 
pas trop bosselée, le premier terme est proportionnel a l’étendue 
de cette surface ; il doit étre négligeable devant le second pourvu 
que la région A ne soit pas trop petite ou n’ait pas la forme d’un 
feuillet. Dans ces conditions, on retiendra que 


(73) An? = — n(N — n)b,o/N 
ou 
(74) An? = n(N — n)[4 + f e— 1) dw,]/N ; 


Ornstein et Zernike avaient trouvé un résultat de cette forme 
dans le cas ot NV est beaucoup plus grand que n, c’est-a-dire quand 
il est permis d’écrire 


n(N —n)/N =n. 


I] nous reste, pour terminer l’étude du fluide homogéne, a 
calculer },,, OU d,, en fonction des données statistiques. Mais 
c’est une erreur de traiter les problémes relatifs au fluide théo- 
rique en supposant qu’il n’y a pas de champ extérieur, ce champ 
est indispensable pour articuler les raisonnements ; son existence 
entraine que le fluide n’est pas homogéne, il est impossible d’avoir 
alors, comme nous l’avons fait ci-dessus, une Opinion empirique 
sur les propriétés de },, aux grandes distances ; nous ne parvien- 
drons 4 élucider la question des fluctuations en densité d’une 
maniére assez générale que par l’application stricte et minutieuse 
des relations fondamentales. 


8. Equilibre thermodynamique. — La répartition des molécules 
d’un fluide placé dans un champ de forces extérieur, invariable 
et maintenu 4 température constante a déja été examinée dans 
la brochure [11] ; les termes en 4 /N avaient paru alors négligeables 
et cette ancienne analyse est insuffisante pour notre objet actuel ; 
une méthode différente va nous permettre de combler cette lacune 
tout en retrouvant les résultats déja connus ; le principe en est 
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du a Ursell [25], il a été repris par Fowler [26] dans sa mécanique 
statistique. Ces auteurs ne paraissent avoir envisagé que les pro- 
blémes relatifs 4 un gaz modérément comprimé non soumis a un 
champ extérieur. L’énergie totale du fluide est 


1 
(75) E = Zz m(uz” + 9:3 + ,*) + 2f, + U2hfie; 


m est la masse d’une molécule, f, le potentiel d’ou dérive le champ 
de forces extérieur ; f, comprend notamment un champ répulsif 
trés localisé mais intense, qui empéche les molécules de traverser 
la paroi ; cette interprétation mathématique d’une paroi physique 
est certainement trés inexacte, mais indispensable, faute de mieux, 
pour que les équations forment un tout complet. f,, est l’énergie 
mutuelle de deux molécules, fonction de la distance r de leurs 
centres ; cette énergie s’annule quand la distance dépasse le rayon 
d’action moléculaire L, trés petit devant les dimensions du réci- 
pient ; elle devient infinie, et positive, Jorsque la distance r est 
inférieure au diameétre moléculaire c. Dans (75) la premiere somme 
est l’énergie cinétique Ec, elle comprend JN termes ; le reste est 
Pénergie potentielle Ep, la seconde somme comprend NV termes 
et la derniére en comprend autant qu’il y a de couples de molé- 
cules, ¢’est-a-dire : N(N — 1). 

Eliminons les vitesses des formules qui nous seront utiles ; 
définissons d’abord ¥, par la relation 


(76) [exp (= EaMkT Aypa08 = 13 
introduisons alors 
v, = ares 

on a 

J exp. (vp — E,)/kKT]dy2...5 = 1 
ou 
(77) exp (— ¥,/kT) = if exp (— E,/kT) dey2.w 3 
considérons 


(78) Dp = exp [(% — E,)/kT] 
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nous avons 


(79) y= Wf Dy dees 


(80) Yo = N(N — 1) | rds. 


ce qui simplifie les formules (7) et (54). 
Il est commode, pour alléger I’écriture, de poser 


(81) 8&1 = exp(—fi/kT) agp. = exp (— f,a/kT) — 1; 


a Vintérieur du récipient, si le champ extérieur est faible, g, 
différe peu de 1; a l’extérieur, au contraire, g, est nul ; quant a 
812, c’est une fonction de r qui est nulle pourr> L,qui est égale 
a — 1 pour r<o et qui reste finie pour les autres valeurs de r 
Au total, cette fonction n’est différente de O que pour un trés 
petit domaine de valeurs de r. Nous pouvons écrire : 


(82) exp.(— E,/kT) = gige--- gy- (4 + S12)(1 + gis)- ++ (4 + goe_ay); 


Dans la suite, le produit des parenthéses sera désigné par Ajo3...y- 
Développons ce produit, puis ordonnons les termes d’aprés le 
nombre de points, M,, M,, My, ..., qwils font intervenir, en com- 
mMengant par ceux qui en comprennent le moins. II vient 


1 1 1 
(83) Aye.y = 41 + 5 =2gi2 + Ezx( 5 812813 + é E126:s800) 4.ses 


dans les sommes multiples, toutes les combinaisons des indices 
sont 4 retenir et il faut distinguer par exemple 812 Ct go, 5 avec ces 
conventions, la somme double comporte N(N — 1) termes, les 
sommes triples en comportent N(N — 1)(N — 2) et ainsi de suite. 
On pose encore : 


(84) Vy = [ gradu, Ng, = iV; 


d’aprés les remarques faites sur la fonction g,, dans le cas d’un 
champ de forces peu intense Vp se confond presque avec le vo- 
lume V et j,, si le fluide est quasi homogéne, avec la densité ve 
Utilisons (82) et (83) pour transformer (77) ; on trouve que 


: 1 
(85) exp (— $,/kT) = Vy" (4 es (1 = a suites 


you 2\ (4 1 a 
=f (1 — x) (1! me x) f (5 812813 + § €is8:s80a) sjadsthoya +++]; 
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le crochet sera désigné par A(J) ; il se présente comme une série, 
limitée, de termes trés nombreux ; comme pour tous les déve- 
loppements analogues de la théorie cinétique, on pourra ne con- 
server que les premiers termes, parce que les termes sont de degré 
croissant, par rapport a j, c’est-a-dire, 4 peu prés, par rapport a 
la densité et parce que les facteurs qui figurent dans les intégrales 


&12 812813 812813823 


deviennent nuls dés que deux des centres de molécules intéressés 
s’écartent l’un de l’autre. 


Essayons de transformer v, donné par (79) comme nous venons 
de faire pour (77). Nous devons ordonner différemment les termes 
de Aj,.... ; considérons un terme quelconque obtenu en effectuant 
les produits de parenthéses, par exemple 

8128138348 67868889 > 


nous le partageons en deux facteurs, l’un qui s’annule nécessai- 
rement quand l’une quelconque des molécules qui y figure s’éloigne 
suffisamment de M,, l'autre qui au contraire ne s’annule pas né- 
cessairement ; ici les deux facteurs sont respectivement 


812813834 et 867868889 > 


le premier facteur sera appelé facteur lié 4 la molécule 1 ; rappro- 
chons tous les termes qui contiennent le méme facteur lié, leur 
somme est 


(86) £12813834 ° As5e67---n 3 


Nous ordonnons ensuite les termes ainsi groupés d’aprés le 
nombre de molécules qui figurent dans le terme lié : 


1 
| Aya... = Agg..w + Dgi2Asa.w + =2( 5 £12213 + 812823 


1 1 
+ 5 Ereerages ) Aas + 2sx( giageaga + 5 812823824834 + 5 8128138 24834 


1 1 
(87) a : 2192132148 24834 + 5 2122138238 24834 + 6 2122139148238 24834 


1 1 1 
+ £12213g23224 + G 612813814 + 5 812813814823 + g12gi3g24 + j guegesees 
Apes es 


la somme simple comporte N —1 termes, la somme double 
(N — 1)(N — 2), la somme triple (NV — 1)(N — 2)(N — 3) et 
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ainsi de suite ; les indices 2, 3, 4,... N peuvent permuter de toutes 
les maniéres possibles. (79) combiné avec (78) et (85) nous donne 


alors 
/ 


1 : 
4 = Ala —1) + AW —2(1—F) [ sraieden 


(88) + A(N — 3)(1 — “ale — Hf (5 sae + £12823 


1 ae: 
Oo 5 EreBisBe | JeJsd%e3 + +++] /A(N). 


Négligeons désormais les termes en 1/N? devant les termes en 
1/N ; il est alors légitime de poser 


4 . 
an S(S—1)(9 —4)] 


S est un entier petit devant V, § un nombre constant dans le 
fluide ; posons encore 


(89) q = Ji:A(N — 1)/A(N), 


il vient 


1 
(90) 4 =i + (1 mer wv) f eretede 


(88 bis) A(N—S)/A(N) = [A(N —1)/A(N)]*[1 — 


Ghd 1 1 
a (1 — y°) ip (5 812813 1 81283 + 5 B1281082s fata +--+]. 


Malgré les remarques qui suivent Péquation (85), cette série 
d’intégrales (dont le nombre de termes est fini) est trés peu con- 
vergente. Pour le montrer supposons que le champ extérieur, en 
dehors du champ dé a la paroi, se réduise au seul champ de pe- 
santeur ; négligeons les termes 4 /N. Choisissons par exemple le 
point M, 4 1 millimétre de la paroi ; essayons de limiter la série 
aux premiers termes ; qj, go, g3, ... ne dépendant que de l’alti- 
tude des points M,, Me. Mae respectivement ; il en résulte 
que les intégrales comme 


if 8129219 i 8128139293dW95, 


ne dépendent de méme que de l’altitude de M,. Ainsi nous pré- 
voyons que la densité v n’est fonction que de laltitude dans un 
fluide soumis a la pesanteur ; ce n’est pourtant expérimentalement 
pas exact a 1 millimétre de la paroi a cause des phénoménes ca- 
pillaires. Au moins pour traiter de ces phénoménes nous serons 
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obligés de prendre dans la série des termes de rang trés élevé de 
maniére quils fassent intervenir le champ au voisinage de la 
paroi:al’intégrale de rang p correspond un domaine d’intégration 
qui s’étend jusqu’a une distance pL du point M, ; il faudra done 
tenir compte de termes de rang p tels que 

pL = Imm. 

Malgré tout, les termes de rang comparable a N sont absolu- 
ment négligeables. 

Le résultat (90) est donc peu maniable. Nous allons le mettre 
sous une forme bien meilleure de la maniére suivante : |’équation 
(90) peut étre considérée comme une équation intégrale ot 
q(x, y, z) est une fonction inconnue ; on résout donc cette équation 
par itération en prenant comme premiére approximation 


y=) 
on néglige les termes en 1/N. En deuxiéme approximation il 
vient 


(91) qa = [1 —{ £12%2d4], 


et finalement, en quatriéme approximation, 


4 4 
{— al £12%2dw_ — 4p (— 5, 812813 = 5 BoSiafan )¥2¥sden 


1 | 4 
= f (5 £r2813824834 + 5) 2102138 14824834 + 5) 8128138238 24834 


ots ee 





(92) ‘ 
ae 6 2198138 14823824834 + 6 612813814 


4 
re, BraSeaSarfas) 2¥0"sd a0 sk ae 


. 


\ 
On obtient, semble-t-il, une relation un peu plus simple en pre- 
nant le logarithme des deux membres. En rapprochant (24) et 
(81) on a, en désignant par K une constante, 
(93) log q = K — f|kT, 
et 





/ 63 4 
K — f,/kT = log “4 —| £19Vgd =| 5 212813823%2"34 23 
| 1 4 
(94) =|. 2108138148238 24834 1 5 6128 13814224834 + 5 812813823824834 


4 
3s 5 Saf is€kas) 09" a0 
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Je pense qu’au contraire de (90) la série du second membre 
de (94) est rapidement convergente ; d’abord elle est plus simple, 
chaque intégrale comporte sensiblement moins de termes ; 
d’autre part, et c’est 14 un renseignement bien indirect et peu ma- 
thématique, elle permet, avec peu de termes, de rendre compte 
des phénoménes capillaires, en particulier de l’existence de couches 
capillaires non planes dans un champ & surfaces de niveau hori- 
zontales. 

Je n’insiste pas sur ce point parce que la formule (94) est déja 
un peu en dehors de mon sujet. Je rappelle seulement comment 
Péquation d’état se déduit de (94) : on suppose d’abord que sous 
action du champ extérieur, dont les surfaces équipotentielles 
sont des plans horizontaux, le fluide reste homogéne a I’échelle 
moléculaire ; dans les intégrales on peut confondre Ya, Vg, +» AVEC 
vy, On a done : 


1 : 
K — f,/kT = logy, — vy J sudo, eet) 128 f esa8iagerdon Soa ae 


ou, en dérivant par rapport a z et en supprimant l’indice 4 14 ov il 
devient inutile, 


— (af/o2)/kT = ov/az(1/v —f sudo, a +f ersfisteos hays 


mais la pression au niveau z est donnée par la formule 


p= | v.0//32.dz ; 
4 


il vient done en intégrant 


1 1 
(95) P/kT = vy — 5 | edo, = 5 | Satiadeadn SET Ess 


C’est l’équation d’état découverte par Ursell [25] [26] et quia 
fait Pobjet d’une étude récente de Mayer [79]. 

Par des procédés analogues 4 ceux qui ont servi pour y,, calcu- 
lons maintenant ¥42 d’apreés la relation (80). Pour éliminer ‘ies for- 
mules inutiles, nous donnerons de la densité simultanée deux éva- 
luations distinctes, la premiére obtenue en négligeant les termes 
en 1/N, valable aux petites distances r ; la seconde qui tiendra 
compte des termes en 1/N, valable seulement aux nie dis- 
tances. Autrement dit nous allons calculer : 19 a,, 20 6. quand la 
distance r est grande (voir (67)). Ps = 
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Pour le premier cas, on écrira comme en (87) 

Ajes3...5 = Aga...x + 2(2:5 + £o3 + £13803) A 45... 

+ =z (5ernlse 18 : 813814834 + 813834 + ; 823824 + ; 823824834 1+ 8038 34 

+ 813824 + 813823834 + 213814823 + 819823804 + 813824834 + 813814824234 

+ 813823824834 + 5 BraBraBesBo4 3% : BrabsaBasBes8e ) Ase. 4] ; 





(96) 


les sommes sont ordonnées d’aprés le nombre d’indices permutables 
qu’elles font intervenir ; dans la somme simple, c’est l’indice 3 
qui peut étre remplacé par n’importe lequel des indices It rat is 
dans la somme double, ce sont les indices 3 et 4 qui peuvent per- 
muter entre eux et avec n’importe lequel des indices Dy Ose IN 
A la somme simple correspondent, avec l’approximation étudiée, 
N permutations ; 4 la somme double N?, et ainsi de suite .(80) 
donne alors 


ie = 9172 ©XP ( — fi2/kT) 4 +f (es + 823 + 813823)93du, 
(97) ; / 
| 13 if (5 813814 + 3 Bia8re8o1+- Matadon ro ; 
on remplace q par son expression tirée de (92), on ordonne les inté- 


grales obtenues d’aprés le nombre de points M qu’elles font inter- 
venir et on obtient le premier résultat cherché 


Vyg=%y%2 EXP (—fj2/kT) A Hi £13823¥3d3++ { (2138 248sat 19814824838 





(98) 
4 4 
| + 812823824834 + 5 813814823824 + 5 8138148228248 04)¥2¥sdgg + - - | 


Si les points M, et M, sont éloignés, tous les produits comme 


£13823 &is8esS3a 0 
sont nuls, f,, est nul aussi, et (98) se réduit a 
Mt Vig bt 

Quand M, et M, sont, au contraire, voisins, les points M,, 
M,,... ne peuvent d’ailleurs s’en écarter sans que les mémes 
produits ne s’annulent si bien que les intégrales ne font intervenir 
que le voisinage immédiat des points M, et M,. On doit considérer 
la série d’intégrales comme assez rapidement convergente. 

Les formules (94) et (98) et autres analogues sont assez difficiles 


a développer quand on fait choix d’une loi de forces particuliére, 


Yvon I. 3 
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mais elles ont le mérite de la rigueur ; elles s’appliquent aux mi- 
lieux gazeux, liquides, 4 la zone de séparation entre un liquide 
et le gaz qui le surmonte. II n’y a aucune raison pour qu’elles ne 
s’appliquent pas a l'état cristallisé. Jusqu’a présent, elles n’ont 
qu’un intérét théorique. Nous allons maintenant arriver avec le 
calcul de b,, 4 une relation immédiatement utilisable et valable 


quelle que soit la loi de forces. 

Reprenons l’identité (96) en Supposant que M, et M, sont éloi- 
gnés l’un de l’autre d’une distance grande devant ZL. Alors, des: 
termes tels que 

813823 81382383488 
sont nuls ; (96) se simplifie 


4 1 
99 Ayjq3-+-x = Agg..n + = (213803) A g5-+-x a 25 813814 a 5 813814834 + £isB 34 
(99) 
1 
+ 5 823824 5 823824834 + Sas83q + S180) Asem + ---]; 


il est commode de poser 


1 1 
8134 = 5 &13814 + 5 813814834 + £13834 
4 + 
(100) 8234 = 5 S23824 ct 5 S238248s4 + Se38s4 


1 
81345 = 813834845 1 5 81s8sa8s5845 +° °° 
(99) devient alors 


Aj93...4 = Aga...x a 2233A 45... a = £03Aas..-n a 22 213804A56..-x =. abel: 
+ XE... 2 813---pfa(p+1)---sA(s41)(94.2)-09 i: 


la somme multiple d’ordre (P, S) comporte 
(NV — 2)(N — 3)-..(W— 5 —1) 
termes équivalents, c’est-a-dire, avec une approximation suffi- 


sante en admettant que les seuls termes utilisés correspondent 4 
une valeur de S petite devant N, 


NU — : S(S — 1)/N}. 
(80) donne alors 


/ 1 3 
“12 = 92 ee n° i (1 i n°) f (es na 823)93dw, tema 
1 
(101) + [1 — 565(S — 1)/N] { 813++-P9374---Godwg4...p X 


{icon eon eete of 9sdop 4.4) (P4.9).008 = . 


35 FLUCTUATIONS EN DENSITE I-35 


Ici une remarque un peu délicate : puisqu’on néglige des termes 
ou une chaine continue de molécules, chacune d’elles se trouvant 
dans le rayon d’action d’au moins une des autres, relie la molé- 
cule 1 4 Ja molécule 2, a-t-on le droit d’appliquer le résultat (101) 
lorsque M, et M, sont distants par exemple de 1 millimétre ? Un 
droit analogue n’existe pas pour l’équation (90) ; mais il faut 
faire attention que la chaine de molécuies considérée a la suite 
de (90) a une extrémité libre soumise 4 la seule condition de longer 
la paroi ; ici la chaine de molécules a une extrémité fixe en M, 
les domaines d’intégrations relatifs 4 cette derniére chaine sont 
donc beaucoup plus restreints maintenant et je crois pouvoir 
admettre que la solution (101) est suffisante. 

Récrivons |’expression (90) de v 4 l’aide des notations (100) ; 
il vient 


1 
\" = q L + (4 — 7°) [ sugedos + 


(102) : 
-- (4 —— 5 P(P wo 10) [ Bie-otte = ‘Jpdwo...p +: rity | ~ 


formons le produit’ v,v, 4 l’aide de cette derniére expression, on 
écrira pour plus de clarté 


4 
juan t+ (13) f sistadon + | 


1 
., == Ys 1+ (19) [ seed oe € 


une fois le produit effectué, on peut permuter les indices de ma- 
niére que le terme général de v,v, soit le méme que celui de vy, 
au facteur entre crochet prés ; dans le cas de v,v, ce facteur doit 
étre remplacé par 
1 
1 — 5 0P(P —1)N — 598(S — 1)IN ; 
la différence v,, — vv, a alors un terme général du méme genre ; 
le facteur entre crochet doit étre remplacé par 
— 6P(S — P)/N; 


cette différence a donc la forme d’un produit ; si on pose, 


(103) vb) = 902 | aes 2 [erated ru 3 { eratstedos aes 4, 
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on a, a l’approximation cherchée, 


(104) Yy2 — Yy%2 = %y¥g0yo/N = — ¥4%Q04b5/N 
ou 
(405) Dis — bbe, 


formules valables seulement si ry. n’est pas petit. b(M,) apparait 
comme une expression qui ne dépend, une fois le fluide donné, 
que des coordonnées du point M,. 

Ornstein et Zernike [27] ont cherché a résoudre des problémes 
analogues 4 ceux qui viennent d’étre traités ; essentiellement, ils 
se sont efforcés de calculer une fonction analogue a Yy2 dans le cas 
d’un fluide homogéne. Leur résultat est fort différent de (97) 
quoique le point de départ soit le méme. II] me semble que certains 
de leurs raisonnements manquent de rigueur : par exemple ils 
découpent le fluide en petits éléments de volume trés voisins, 
contenant seulement quelques molécules ; ils admettent ensuite 
que les forces exercées par les molécules d’un élément sur les 
molécules de l’autre se résument en une force unique dirigée 
suivant la droite qui joint les centres des deux éléments : cette 
affirmation est certainement inexacte quand les deux éléments 
se touchent. Un exposé clair et concis des idées d’Ornstein et 
Zernike a été donné par Didlaukies [28]. 


9. Nouvelles formules de fluctuations. — La densité vy, est fonc- 
tion : 1° de M, 2° de la température, 3° du champ extérieur, 
4° de N. Calculons d’apzés (102) 


ov, / oN ? 


en négligeant les termes en 1/N ; de (88 bis) et (89) on déduit 
aprés quelques calculs 


2 log g,/aV = 6/N; 


et on voit immédiatement que 


ae 
(106) vb, = NO 2 av /aN; 


dv,/9N est le taux de variation de la densité moléculaire au point 
N, quand on ajoute une molécule dans le récipient. Comme 


Jae, = N, 
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ona 

(i Nev, /aN dw, = N 
et 
(407) if ¥,5,dw, = NO73. 


La relation (70) que je récris, 


J bude =—1 — | les — 1)dw, 


(408) et nee eee fh Vel dig — ddtng: 


devient donc 


en commettant un erreur insignifiante, parce que l’on confond 
b,, avec — b,b, méme pour r,, trés petit. 

Les raisonnements qui ont conduit de (66) a (73) peuvent étre 
repris sans modification maintenant dans le cas d’un fluide 
inhomogéne et donnent pour le carré moyen des fluctuations 
lexpression 7 


oO 


(109) An? = 073 { ¥,5, do, - | v4, deo,/ if v,b,dw, 
V—A vy 


vA 


résultat valable seulement dans le cas ot le volume A contient 
en moyenne un grand nombre de molécules, ot la surface qui li- 
mite ce volume est 4 peu prés lisse. Si le volume V est beaucoup 
plus grand que le volume A, on a sensiblement 


(110) An? = Ff hide, 
A 


10. Fluctuations dans le champ de pesanteur. — La fonc- 
tion b(M) définie par (103) peut s’exprimer trés simplement 
quand le champ de forces extérieur est uniforme, ce qui est le cas 
du champ de pesanteur défini par le potentiel 


{(M) = mgz; 


m est la masse d’une molécule, g est l’accélération de la pesanteur 
en valeur absolue, z est l’altitude. 
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Calculons dv/dz d’aprés l’équation (90) ; au préalable on néglige 
les termes en 1/N. On remarque par exemple que 


(111) = f sutdes 2] £19(249/022) dw, 


et de plus que 
kT - aq/az = — mgq; 


il vient 
kT - ay,/02, = -— mgljy + dir f eshte ae Bir { Suoinisden +-+-];5 


eae yee 
on reconnait dans le crochet l’expression (103) de v,b,9 2. Si bien 
que nous arrivons finalement 4 la relation 


is 
(412) kT - dv,/02, = — mgy,b,9 2 
qui donne b sous une forme trés maniable. 

Appliquons ce résultat au calcul du carré moyen (110) ; pour 
simplifier, choisissons comme volume A un cylindre vertical li- 
mité par deux faces horizontales de section S ; supposons de plus 
que les couches d’égales densités soient horizontales a l’intérieur 
du cylindre, alors v est fonction seulement de z. On a 


An? = — (SkT|mg) f ov/osae 


Pintégrale est prise sur la hauteur du cylindre ; désignons par v, 
la densité a la base, par v, la densité au sommet, on a 
(113) Ant = SkT 24. 
mg 

Introduisons dans cette formule des grandeurs plus familiéres : 
remplacons k par R/N,, v par No, m par M/N,; M est la masse 
moléculaire, p la densité en nombre de molécules gramme par 
unité de volume. II vient 





Ant — od becowes y| 
(114) An? = NRT Me S. 


La formule (114) est la réponse au probléme que nous avions 
posé a la fin du § 5. Elle est satisfaisante parce que: 1° elle donne 
pour le carré moyen des fluctuations une valeur finie méme au 
voisinage du domaine critique, 2° elle satisfait a la relation (43) 
comme d’ailleurs la formule préliminaire (110), 3° elle redonne, 
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lorsque le fluide est loin de ]’état critique, le résultat de Smolu- 
chowski (40). En effet, loin de l’état critique, en supposant que 
dans le cylindre A le fluide ne présente qu’une seule phase, li- 
quide ou gazeuse, on a sensiblement d’aprés (47) et (48) 


Pe — Py = — (aP/2p)—?- Mgo(z_ — 2) = — Mgr*S(z, — z,); 


dans ces équations ¢ représente la densité moyenne. II vient 
d’aprés (114) 


An? = NoRT¢*25(2 — 22), 
mais 
n = NopS(2, — 2) 
est le nombre de molécules contenues dans le cylindre, si bien que 
Pon retrouve la formule 
An? = RT Gon. 


Imaginons maintenant que la pesanteur tende vers zéro. La 
formule (114) perd toute signification a la limite, si la température 
est la température critique, parce qu’alors Ja densité n’est plus 
une fonction de l’altitude z. Les variations de la densité 4 l’inté- 
rieur du fluide ne sont plus commandées seulement par le champ 
de pesanteur, mais encore par le champ 4 la paroi, le phénoméne 
est sous la dépendance des forces capillaires introduites par Ro- 
card [24]. Pour éliminer les parois, on peut imaginer qu’elles 
s’éloignent indéfiniment en méme temps que le champ tend vers 
zero ; ’équation du probléme fait intervenir l’expression 


07y/oa2 + 92/ay? + 0a°v/02", 


il en résulte que l’effet de la paroi ne diminue pas quand elle 
s’éloigne ; il est impossible de la négliger. La formule de fluctua- 
tion que nous pouvons retenir alors est (109) ou (110) : elle montre 
_que les fluctuations restent de toute maniére finies. 


44. Formule d’Ornstein et Zernike. — Dans le cas d’un 
fluide quasi-homogéne, la relation (112) devient d’aprés (48), en 
supprimant les indices devenus inutiles 


wt! 
kTN,Mgr?2 = mgvbs ? 
ou 


1 
(445) b6 ?= RT Ze; 
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mais (108) s’écrit, en supprimant 1a aussi des indices inutiles, et 
puisque le fluide est par surcroit isotrope, 


ee 4 + anv f far — 1)r*dr. 
On a done 
(116) tor f [a(r) — 1]r*dr = RT Bo — 1. 


Cette importante formule est due, avec des notations différentes, 
a Ornstein et Zernike (voir [28]). Ces auteurs l’établirent beaucoup 
plus facilement qu’ici, ils ne considéraient que le cas d’un fluide 
homogéne et n’eurent qu’a combiner (40) et (74). J’en donnerai 
une démonstration directe et qui fait appel au théoréme statis- 
tique de Lionville & propos de la diffusion de la lumiére. Je rappelle 
que a, a été introduit par (67); 4zv a(r)r.dr représente le nombre 
de molécules qui se trouve en moyenne dans la couche sphérique 
de rayon r, d’épaisseur dr, centrée sur une molécule. a(r) est nul 
pour r<o et est égal & 1 quand r dépasse largement le rayon d’ac- 
tion L des forces intermoléculaires ; les valeurs intermédiaires 
de a(r) peuvent étre déduites de la diffusion des rayons X par le 
liquide suivant le procédé donné par Zernike et Prins [15] (la 
méthode de calcul a été précisée récemment par Randall [87)}). 

Dans un gaz parfait, on a 


RTfSo—1=—0; 


pour un liquide la méme quantité est négative ; prenons, pour 
fixer les idées, le benzéne, quoique ses molécules ne puissent guére 
passer pour sphériques. On a alors, en c. g. s., 


R = 83.107 T = 290° 2 = 88.10—12 e = 0,011 
et 
RT6o —4 — — 0,98; 


au voisinage du point critique cette quantité est au contraire 
positive : il semble qu’alors les molécules aient tendance a se 
grouper en amas laissant entre eux un vide relatif. 

La formule (116) doit rendre service dans linterprétation des 
diagrammes de rayons X dus a la diffusion de ces rayons par les 
fluides ; il est curieux quelle ne figure pas dans le mémoire [15] 
de Zernike et Prins sur la question, toutefois ces auteurs y font 
implicitement allusion au cours de remarques qui terminent leur 
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mémoire. Cette omission a entrainé une erreur de Gingrich et 
Warren [29] qui ont cru pouvoir poser 


key i [a(r) — 4]"dr = —4, 


c’est-a-dire sous une autre forme 


1 
of Meu — i)du, = —1; 


alors qu’une relation exacte voisine est 


= | G42 = Ade, = — i, 


Gingrich et Warren ont ainsi utilisé une relation presque exacte 
pour un liquide, mais trés inexacte au voisinage de point critique. 
Il faut reconnaitre que la distinction entre v,, et v,a,, est un peu 
subtile. 


CHAPITRE II 


FLUCTUATIONS EN DENSITE ET EN CONCENTRATION 
DANS LES MELANGES BINAIRES 


1. Définitions. — Deux corps purs sont mélangés dans un réci- 
pient de volume V et forment une ou plusieurs phases ; il n’y a 
aucune réaction chimique entre eux. Le récipient contient NV mo- 
lécules du premier constituant et NV’ du second. 

Une fraction A du volume total contient un nombre variable 
de molécules de chaque constituant. Nous désignons respective- 
ment par n et par n’, les valeurs moyennes de chacun de ces 
nombres et par 


nit \An > Cet® wen Ane 
les valeurs variables. An et An’, sont des fluctuations qui restent 


toujours petites devant n et n’; les valeurs moyennes des fluc- 
tuations sont nulles : 


(147) An = 0 An' = 0, 
I] n’en est pas de méme des carrés moyens 
(418) An? et —_ An®, 
et aussi du produit moyen 
(119) AnAn’. 
Nous nous proposons d’étudier ces moyennes, qui se rattachent 
au probleme de diffusion de la lumiére par les mélanges binaires ; 


nous séparons d’ailleurs complétement cette étude des recherches 
qui intéressent directement la diffusion. 

Nous supposerons d’abord que le fluide contenu dans le volume 
A est homogéne. Les densités moléculaires moyennes, respecti- 
vement v et v’ y sont constantes et ont pour expression : 


v= nA v= n'lA; 
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nous définissons aussi la densité globale en nombre de molécules 
par unité de volume : 


(n + An +n’ + An')/N,A, 
et la concentration 
(n' +- An')/(n + An); 


Ja densité globale et la concentration ont pour valeur moyenne 
respectivement 


(120) p= (n + n')/NoA = (v + v’)/No, 
et 
(121) $=.n'in = v'/vs 


aux fluctuations An et An’ correspondent des fluctuations de la 
densité et de la concentration qui sont respectivement 


(122) Ap/p = (An +An’)/(n + 1’) 
et 
(123) As/s = An'/n' — An/n; 


naturellement on traite toutes ces fluctuations comme des infi- 
niments petits. Les moyennes 


(424) As? Ac? = Apds, 


peuvent s’exprimer linéairement 4 l’aide des moyennes (118) et 
(119), et réciproquement. I] en est de méme pour toute combinai- 
son linéaire des As et Ao. 

Einstein a calculé As* dans son mémoire [17]; il négligeait 
les fluctuations en densité. Raman et Ramanathan [30! ont ensuite 
repris la question en tenant compte de ces fluctuations, leur ré- 
sultat pour As? différe un peu de celui d’Einstein. Les raisonne- 
ments de Raman et Ramanathan, qui ont été résumés par Rocard 
[16], utilisent pour le calcul du potentiel thermodynamique dans 
Pétat fluctué la méthode des parois semi-perméables, qui n’est 
pas trés claire. Je vais reprendre la question en utilisant unique- 
ment des différentiations sires : de Raman et Ramanathan 
les résultats doivent étre modifiés. 

Dans un mélange, la probabilité des fluctuations est donnée 
par une formule qui généralise (26) : 


(125) C exp. [— An?.a?¥/an® — 2An- An’. a?W/anan’ — An’?-2?W/on?)/2kT, 
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C est une constante, VY est le potentiel thermodynamique géné- 
rique ; on a supposé pour simplifier V beaucoup plus grand que A. 
On sait que ¥ est de la forme 


(126) Ws AG(T Hvev ss 
ona 
AeW/an® = a*b/av? —-_Aa®W/anan' = 9° /avev' 


(127) Adan? = a?&/avav’; 


nous désignerons pour abréger ces trois dérivées secondes de ® 
respectivement par G, H, J. La probabilité (127) peut donc s’écrire 


(128) C exp. [— GAn* — 2HAnAn’ — IAn'*]/2 AkT. 


Nous allons exprimer les coefficients G, H, J, 4 aide de para- 
métres ayant une signification concréte. Remarquons l’abord que 
la pression uniforme P qui régne dans le volume A est donnée par 


(130) P =—oW/aA = — © + v-ab/av + v'-aP/av’; 


introduisons les pressions partielles des deux constituants, res- 
pectivement p et p’, définies par les égalités 


(131) f+ kT log p=oav/w f' + kT log p! = a®/av’; 


si le mélange est en équilibre 4 travers une paroi semi-perméable 
avec l’un de ses constituants, si celui-ci est gazeux et se comporte 
comme un gaz parfait, la pression de ce gaz se confond avec la 
pression partielle correspondante. f et f’, sont des fonctions de la 
temperature seulement, qui dépendent du choix de la fonction 
®, toujours un peu arbitraire en thermodynamique ; 
P, p', e sont des fonctions de v, v’, T, mais aussi bien de UE teed § 
c’est ce que nous supposerons ci-dessous. 

Différentions membre 4 membre les relations (130) et (131) en 
laissant le volume et la température constants, il vient 


(132) dP = Gvdv + Alv'dv + vdv’) + Iv'dv! 
op Q 

(133) xr| 22 EEN 5 AP | /p Si Gde sae 
op’ op’ : 

(134) kT [22 as + oF, oP | [p' = Hav + Id’; 


en éliminant dv et dy, entre ces relations on trouve facilement la 
relation classique 


(135) SP Dee 
P08) 2D eae oes 
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et de plus que 


yop | v ap’) _ 
(136) xr enitas 5) 1. 


Mais introduisons maintenant le coefficient de compressibilité 
du mélange 
2 =0 log p/aP, 


et aussi le coefficient de variation de la densité avec la concen- 
tration 


(137) L = log o/as; 

ona 

(138) dy'/v’ — dv/v = ds/s, 

(139) d log p = (dv + dv’)/(v + v’) = Ids/s + 2dP; 


ces relations servent a éliminer ds et d Pde (133) et (134) ; on 
trouve alors, d’aprés (133), 





(140) G = — kTs(a log p/as)/v + kT (a log ploP)( . + aE 


y + y! v 
H = — kTs(alog p/as)/v' + kT (a log p/aP) (. ! 1 ae ib 


et d’aprés (134) 


l 
(1441) H = — kTs(a log p'/ds)/v + kT (a log pleP)( + aE 





1 l 
= — kTs(a log p'/as)/v' +- kT(a log p'/aP) (; ra 7 

On égale les deux expressions de H ; compte tenu de (135) et 
(136) il vient 

4 
v+y’ 

On n’utilise pas davantage les relations (140) et (141), mais a 
Vaide de (133) et (134) on forme l’expression 


Gdyv* + 2Hdvdv' + Idv’, 


l 
kTa log p/aP = + : kToa log p'JaP = ——, — =: 


ou on exprime, grace a (138) et (139), les différentielles dv et dv’ 
en fonction des différentielles dP et ds. I] vient 

ds 2 
(142) Gdv* + 2Hdvdv' + Idv® = — kTv'-(a log p/es)- (2) + ((dP)?. 
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Il est donc naturel de considérer la fluctuation de pression AP 
définie par 
(143) BAP = — lAs/s + A loge; 
cette fluctuation n’a pas une signification physique simple comme 
An par exemple, mais c’est un intermédiaire commode pour les 
calculs. A l’aide de (128), il est facile de voir que 


(144) AP? = kT/AS = RT/N,A2 

(145) APAs = 0 

(146) As?/s? = — 1/[Av'(a log p/as)]; 

on en déduit aisément 

(147) N,AAe?/o? = RT’ — PN,/v'(a log p/as) 
(148) NyAdpAs/ps = — IN,|v'(2 log p/as). 


Raman et Ramanathan trouvaient 
Kee ea 
NoAKA? = RTB. 


Le point critique de miscibilité compléte d’un mélange binaire 
est. défini par les conditions 


(149) 2 log p/ds = 0 

2? log p/ds? = 0; 
contrairement aux conclusions des auteurs que je viens de citer, 
les fluctuations (147) et (148) sont infinies au point critique, 
comme d’ailleurs la fluctuation (146). 

Les formules précédentes pourront servir ultérieurement a 
amorcer une théorie de la diffusion dans les mélanges binaires, 
mais il parait actuellement prématuré d’essayer d’améliorer les 
formules de Raman et Ramanathan sur la question. Bornons- 
nous a signaler les études expérimentales de Parthasarathy [83] 
sur lintensité de la diffusion et la dépolarisation, de Gans et 
Stuart [84] et de Krishnan [85] sur la dépolarisation. 


2. Les mélanges binaires dans le champ de pesanteur. — Les 
différentes formules (146), (147), (148) prévoient pour le point 
critique de dissolution des fluctuations infinies, ce qui n’a physi- 
quement aucun sens ; les mélanges binaires présentent donc une 
difficulté analogue a celle rencontrée 4 propos de la formule (40) 
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pour un fluide pur. A mon avis, cette difficulté n’est qu’apparente 
et s’élimine de la méme maniére que dans le cas d’un fluide pur : 
la pesanteur empéche l’existence d’une phase homogeéne qui vé- 
rifierait les relations (149). 

Cette affirmation est loin d’étre aussi évidente que |’affirmation 
analogue relative & un fluide unique ; il semble au premier abord 
que la pesanteur ne peut avoir aucune action sur un mélange bi- 
naire qui se trouve dans le domaine critique, ce qu’on reconnait 
a Popalescence bleue qu’il présente, parce que ce mélange n’est 
pas plus compressible qu’un liquide ordinaire. Mais nous allons 
voir que la pesanteur peut empécher un fluide binaire de rester 
homogéne non seulement s’il est compressible, mais aussi pour 
une autre raison. 

Voila en tout cas comment il faut se représenter les phénomeénes ; 
raisonnons pour fixer les idées sur le mélange nitrobenzéne- 
hexane normal dont la température critique voisine § = 1994. 
On imagine qu’on fait refroidir trés lentement un long tube ver- 
tical rempli du mélange avec une concentration appropriée qui 
n’est pas trés exactement définie et qui correspond & 50 % en 
masse de chaque constituant. La température initiale surpasse 9 
de quelques degrés et le fluide est alors homogéne ; lorsque la 
température du mélange ne dépasse plus § que de quelques milliémes 
de degré, ’homogénéité cesse, le nitrobenzéne, qui est le cons- 
tituant le plus lourd, se concentre un peu vers le fond ; la con- 
centration varie de maniére continue, depuis le bas jusqu’au haut 
du tube ; lorsque la température baisse encore cette variation 
de concentration s’accentue un peu et en méme temps se localise 
dans une région du tube qui mesure quelques millimétres de hau- 
teur : deux phases 4 peu prés homogénes apparaissent, l’une plus 
riche en nitrobenzéne que l’autre et séparées par une couche ov 
la densité varie rapidement et qui présente le phénoméne de 
mirage. Quand la température atteint enfin 6, la zone de mirage 
prend une épaisseur extrémement faible, dont l’ordre de hauteur 
est une dizaine de fois le rayon d’action moléculaire : ainsi est 
apparue entre les deux phases une couche capillaire. 

Si cette description est vraisemblable, elle ne correspond malheu- 
reusement jusqu’a présent 4 aucune observation réelle : personne 
jusqu’A présent ne s’est soucié d’observer le point critique de 
miscibilité compléte de deux liquides avec la précision que Gouy 
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a apportée dans le cas de]’anhydride carbonique pur. Ces recherches 
seraient considérablement génées par la viscosité des liquides 
qui retarde la réalisation de l’équilibre ; une étude continue du 
phénoméne est impossible : 4 chaque température d’observation 
il est nécessaire d’agiter le mélange pour faciliter le passage de 
chaque constituant d’une phase a l’autre. 

Quoiqu’il en soit, entreprenons l’étude statistique d’un mélange 
de deux sortes de molécules. L’étude analogue déja faite sur un 
fluide pur me permettra de sauter des intermédiaires de calcul 
dont la répétition serait fastidieuse. 

L’ensemble des notations est le suivant : 


premier constituant deuxiéme constituant 

Masse des molécules.......... m m' 
Nombre de molécules........ N N' 
Numeérotées—25 Nes cee are 5 Hy Re ns | LTTE SSN 
Masse moléculaire .......... M M’ 
Elément de volume décrit par 

le centre d’une molécule.... dw, dw, 
Energie dans le champ de 

force extérieur «i, cys ee h ht 
Energie mutuelle de deux mo- 

lécules de méme espéce .... hie fin 
Nombre moyen de molécules 

dans un volume donné A ... n n’ 
Fhuictuations: <0 ee An An’ 


fu est Pénergie mutuelle de deux molécules de natures diffé- 
rentes. 

Les densités ne calculent pas une généralisation des formules 
(7) et (54). La densité moléculaire du constituant n° 1 au point 
Z, ¥,% est donnée par 


(150) y= N ff Dévoeresenn dyi2...5' 5 


D est le coefficient de probabilité ; les intégrations sont étendues 
a toutes les valeurs des vitesses et a toutes les positions des centres 
de molécules 4 l’intérieur du récipient. De méme la densité molé- 
culaire du second constituant au point x,y,z, est 


(151) y= Nf Daoy-onsimy dy12...9! 3 
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il y a maintenant a considérer trois densités simultanées doubles : 
(152) Yo = N(N —1) { f Ddvseexr dy12...n' > 
Vy = ww’ || Ddwes...x nim sey’ dy12..-y’ ’ 


Yn = N'(N' = 1) [ f Pavreommy dy 12... 3 


ces différentes densités vérifient les relations intégrales suivantes : 


(153) f ndes = N J nde, ==. JV’ 


or 


(154) J rasdos = (NV aes 1)v4 J Vy 7d = N'v4 


J odes = Nv, i] Wyndwy = (NV — 1)y, 


A laide des densités simultanées on exprime les moyennes de 
fluctuation dans le volume A, sur lequel il n’est fait d’abord aucune 
hypothése restrictive : 


(155) Dt — fof (ra — wreldoydeg 
AJ V—A 


AnAn' = — ip | (vir — va¥) dodo, 
AJ) V—A 


An? as fey: (Vin — vyp/n) dwydwy 
JA —A 


I] parait vraisemblable a priori que les densités simultanées 
peuvent étre représentées par des séries entiéres en 1/N et 1/N’ 
dont on ne conserve que les premiers termes : 


(156) Yao = %%e(@yg + Oyo/N + b'y0/N’) 
Yay = %4¥7(Aay + by;/N + bi;/ ’) 
Yn = y%n(AQn + bin/N + bin/N'); 


les coefficients dy., 1, Qn, Oye, ..., bry dépendent uniquement 
des propriétés du fluide dans le voisinage des points M,, M;, Mg, 
My, ..-- Lorsque la distance M,M; augmente par exemple, ay 
tend rapidement vers 1, 6, et b', tendent chacun vers une limite 
qui n’est pas nulle. Il est d’ailleurs légitime de confondte les b 


avec leurs valeurs limites respectives pour l’usage que nous voulons 


Yvon I. vA 
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en faire. Pour aller plus loin, nous devons nous placer dans le cas 
de l’équilibre thermodynamique avec un champ de forces exte- 
rieur ; on pose 


gi =exp(—fi/kT) gr =exp(—fi/kT) giz = exp (— fax/kT) — 1 
8u = exp (—fir/kT)—1 gin = exp (— fin/kT) — 1; 


on pose aussi 
V,= fede, Vi=[ gde, 
Ng = jiV> N'g1 = JV. 


A(NV) défini a la suite de (85) est remplacé ici par A(N, NV’). 
On pose 

91 = Jr: A(N4,N')/A(N,N’) 

gx = Jr A(V,N'_1)/A(V,N’) 


et on a en négligeant les termes en 1/NV, 1/N’, 
(157) v1 = gill + fe 212q2d2 + of £ugide; 
4 4 
ss af (5 812813 + £12823 + 5 S18isdee )gagsdon 
+ f (Zie811 + Sie8ue + gieger +21u821)q2g:dwe; 
4 1 
+f (Senein + 5 Susi ngin ate gigi )digadorn +-:-] 
LETT | a iS 1 
“ = af cia £ingudey +] 2119101 +f 5 SinSim + gingom 
1 
+ 581 nfimgom)9n9meon om 
+f (gingn + 2m811%1n + finger + £11821)929nden 


4 4 
+f (Senge = 5 8uS21812 > gigs angadon +--+]; 


ces relations sont peu maniables ; on les résout par rapport a 
q, et gq ; plutot que g, et gq: d’ailleurs, on calcule leurs logarithmes 
et on a finalement avec deux nouvelles constantes K et K’ : 


(158) K— fi/kT = log V4 i f £12ved2 — iB £1r1dw; — ; J grzgigesders 
4 ’ 
= { 2122118 21VeV1d. wey mas | Sulinginn1dyy — --- ; 
K'— fi/kT = log oF —| 2m/n1dwy — | girdoy —- 5 f érngemgxmdn 


» 1 . 
ao 7 Ornd11810%1¥1d 4 — 3 { 811812212¥1"2dw32 — +--+, 
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Ces équations, qui nous écartent, pour le moment, des fluctua- 
tions, permettent de discuter la répartition des molécules sous 
Vaction d’un champ extérieur, et nous pensons surtout au champ 
de pesanteur. Nous supposons le récipient assez large pour que les 
parois ne jouent qu’un role accessoire. 

Plagons-nous dans le cas ot le champ extérieur ne crée qu'une 
faible hétérogénéité ; alors, dans (158), en supposant que M, 
et M, sont voisins, on peut confondre vo, vs... avec vy eb Vy, Ym. 
avec v;; les seconds membres de (158) ne dépendent alors de 
M, et M; que par l’intermédiaire de v, et ». 

Posons 


4 
(159) B(v,,¥1) = v,(log vy; —1) + (log »—t)—t f eredo, = wef bud 
4 2 : 1s f ae 
meres Bs £1ndwy — OF J 21281382343 — 91 “1 £128 118191 


t 1 
=< rtf Sugingendosn— 5 | eeneemfamonm — ses} 


On vérifie immédiatement que (158) s’écrit 
K — f,/kT = aBloy 
K' — f, [kT = oB/oy; 
dans un champ de forces verticales les densités ne dépendent que 
de z; on élimine les constantes K en dérivant : 
(160)  —oaf,/oz = (22B/av,2+ avy/dz + 2°B /avd»- dv /dz)kT 
— af, /az = (22B /avyovy + av, /dz + o°B Jove + dvy/dz)kT. 


La pression au niveau z est donnée par la formule 


P= { [v4(af; /22) + ve(ofr/2z) dz 


v2 
ce qui s’écrit d’apres (160) 
(161) P = v,2B/av, + vepB/o, — B; 


en rapprochant ce résultat de (130) on voit que les fonctions 
D(v, v’) et B(v, v’) sont identiques 4 une fonction linéaire de v 
et v’ prés. 

Kirkwood [31] avait déja étudié le probléme du fluide homogéne 
et donné une équation d’état équivalente a (161) explicité d’aprés 
(159). 

Désignons maintenant les densités uniquement par les nota- 


52 FLUCTUATIONS EN DENSITE I-52 


tions v et v’ ; appliquons (160) au cas du champ de pesanteur, 
adoptons les notations qui ont servi pour (120). Il vient 
(162) — mg |kT = G- dv/dz + Hav' /dz 

— m'g/kT = H. dv/dz + I dv’ /dz; 
(133) et (134) nous permettent d’introduire d’autres dérivées et 
d’écrire 


oapds apdP 
oe) ape sh rer] 22 dz + Pal 
. ap’ ds | ap’ =| 
tal rere] as dz + 9P dz|" 


Ces équations montrent comment la pesanteur agit pour faire 
varier la concentration et la pression avec l’altitude. On en déduit 
en particulier 


(164) dlogs /dz = gy[ MI |v’ + M'l/v—(M + M’) |(v + v')]/R2T%(a logp /as); 


lorsque le mélange tend vers l’état critique, dp /ds tend vers 0, 
il en résulte que ds /dz augmente beaucoyp, le fluide finit méme 
par devenir trop peu homogéne pour que l’analyse précédente 
a partir de (159) reste valable : c’est alors qu’apparaitrait la zone 
de mirage. 

Quant 4 la variation de la pression avec Paltitude, elle reste 
finie, ce qui n’a rien de remarquable dans le cas d’un fluide qui 
ne tend pas 4 devenir infiniment compressible. 

Nous avons précisé la structure théorique du fluide dans lequel 
nous voulons étudier les fluctuations ; revenons au calcul de ces 
derniéres, d’abord en supposant simplement que le fluide est en 
équilibre thermodynamique, sans rien spécifier sur les particu- 
larités du champ extérieur. 

Dans ce but, j’ai calculé les trois densités simultanées doubles 
dans le cas des grandes distances. J’ai obtenu trois expressions : 


(165) 9 = vyvg — &482D /W — mange /W + (em9 + ean) L/W 
Yar = ¥y¥y — €y71D /W — en, /W + (8 + nya) E /W 
Min = Wy — §&yC /W — agD /W + (ran + en) E/W; 


ces formules confirment les développements (156) : dans chaque 
second membre le premier terme est, comparé aux suivants, du 
méme ordre de grandeur que NV ou N’ par rapport a 1. Lorsque 
le fluide est donné, C, D, E, W sont des constantes et les quatre 
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fonctions 7, 4, m1, « ne dépendent que du point M, ou du 
point M,. 


Les for mules suivantes précisent le sens de ces nouvelles gran- 
deurs : 


(166) y= +f verlau — 1)dw, or FY + f vora(tn — 1)dwg 


"1 = f vertex — 1)dw, 11 =f vonten — 1)do, 
(167) C= fede D = [ ado, B= [do = = [ nude, 


W = CD — E’; 


les moyennes de fluctuations dans un volume A petit devant le 
volume total mais grand devant les dimensions moléculaires 
sont d’aprés (155) : 


(168) An? = i] edo,  Andn’'= { ndw, An®= if edu. 
A A A 


Ces relations se précisent lorsqu’on suppose que le champ de 
forces extérieur est le*champ de pesanteur. Dans ce cas, que le 
fluide soit homogéne, ‘séparé en deux phases ou qu’il comporte 
une zone de mirage, on établit as relations suivantes : 


(169) g[me, + m'ny] = — kTov,/0z 
g[many + me] = — kToy,/az; 


dans ces équations qui généralisent (112) les points M, et M, 
doivent étre regardés comme occupant la méme position. 

Choisissons comme volume A un cylindre vertical de section S. 
Les densités sont a la base v, et vy, respectivement, 4 la face supé- 
rieure v, et v; ; dans tout le cylindre les surfaces d’égales densités 
sont horizontales ; en combinant (168) et (169) il vient : 


(170) g[mAn? + m’AnAn'] = kTS(v, — v4) 
g[mAnAn’ + m’An®] = kT'S(v_ — 1); 


il ne m’a pas paru possible de trouver une troisiéme relation 
linéaire simple entre les moyennes, qui en aurait permis un calcul 
complet. On peut vérifier que, lorsqu’on suppose le fluide quasi- 
homogéne, les moyennes calculées a l’aide de (144), (145), (146) 
satisfont 4 (170). Mais les relations (170) sont plus générales que 
(144), (145), (146) ; elles restent finies quel que soit ]’état du fluide. 
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Il semble done prouvé, avec une lacune il est vrai, que la théorie 
des fluctuations dans un mélange binaire ne comporte pas de diffi- 
culté au voisinage du point critique de miscibilité pourvu que Von 
tienne compte de la pesanteur. 

Je termine ce chapitre par quelques remarques d’ordre expé- 
rimental. 

D’aprés (164) ou d’aprés les formules plus générales (158) il 
est clair que les phénoménes au voisinage du point critique dé- 
pendent de l’intensité du champ de forces extérieur ; or, il est 
facile de réaliser un champ d’inertie beaucoup plus intense que 
le champ de pesanteur en soumettant le mélange, placé dans des 
tubes appropriés, a l’action d’une centrifugeuse. Pendant le 
fonctionnement il tend a se produire, d’un bout 4 l’autre du tube, 
une variation de concentration plus importante que sous le seul 
effet de la pesanteur ; en observant le liquide immédiatement 
aprés arrét de la centrifugeuse, la variation de la concentration 
n’a pas le temps de diminuer et elle serait facile 4 mettre en évi- 
dence si elle était assez forte pour provoquer la séparation du 
liquide en deux couches. I] est malheureusement impossible de 
prévoir si cet effet existe dans un large intervalle de température 
ou seulement au contraire dans un intervalle étroit. J’ai fait un 
essai avec le mélange nitrobenzéne-hexane commercial et une 
petite centrifugeuse; il n’y a eu aucun résultat, ce qui ne contredit 
guére la possibilité du phénoméne envisagé parce que la tempé- 
rature varie beaucoup pendant le fonctionnement, de prés de 41°. 


Une autre maniére de modifier les propriétés du milieu consiste 
a le soumettre a des vibrations sonores ou méme ultrasonores ; 
on pouvait s’attendre, dans le mélange opalescent, 4 des varia- 
tions locales de concentration sous l’effet de l’onde élastique plus 
importantes que dans un mélange ordinaire (phénoméne envisagé 
depuis longtemps dans le cas des gaz par M. M. Brillouin, voir 
[32]) et susceptibles d’avoir un effet sur la vitesse de propagation. 

Lucas et Biquard d’une part, Debye et Sears de l’autre, ont 
découvert (voir [33]) le phénoméne de la diffraction de la lumiére 
par les liquides parcourus par les ultra-sons qui permet de mesurer 
aisément la vitesse des ultrasons dans ces milieux. M. R. Lucas 
a bien voulu, et je l’en remercie ici, entreprendre de mesurer cette 
vitesse dans un mélange opalescent [89]. Les ultrasons, émis par 
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un quartz de fréquence §.10°, sont produits dans une cuve pleine 
de toluene. Un tube en pyrex muni de 4 fenétres planes en verre 
assez mince et paralléles deux 4 deux, contient le mélange opa- 
lescent ; il est plongé dans le toluéne ; par une fenétre les ultra- 
sons parviennent dans le mélange qui est traversé par un faisceau 
lumineux perpendiculaire au faisceau sonore. La lumiére émergente 
est observée visuellement dans une lunette dont loculaire est 
pourvu d’un micrométre. Le mélange étudié, essence de térében- 
thine-aniline, dont la température critique était voisine de 19°2, 
a été examiné a des températures croissantes obtenues par simple 
réchauffement du bain de toluéne au contact de la salle ; la tempé- 
rature était lue directement dans le mélange avec un thermométre 
au 1/10. La température était sans doute définie a 0°05 prés. 
Les observations ‘furent faites en dessous de la température cri- 
tique, sur le fluide séparé en deux couches, puis sur le fluide homo- 
géne. Sur un intervalle de 2 a 3°, la vitesse des ultrasons n’a pas 
varié, dans la limite de la précision des mesures, 10 % environ. 
Il n’y a donc de ce cdté aucune anomalie. Par contre, le mélange 
et aussi le mélange hexane commercial nitrobenzéne, a paru 
trés absorbant pour les ultrasons: c’est 1a une mesure qui nécessite 
un montage plus complexe et qui reste a faire. 


CHAPITRE III 


REPRESENTATION ONDULATOIRE 
DE L’AGITATION THERMIQUE 


J’ai rappelé au § 4 la méthode inaugurée par Einstein pour re- 
présenter l’agitation thermique a l’aide d’une série de Fourier. 
Je vais montrer maintenant comment il est possible d’arriver 4 une 
représentation de ce genre, tout a fait rigoureuse dans le cadre 
de la mécanique statistique classique,ou une déduction mathéma- 
thique satisfaisante prendra la place des hypothéses empiriques 
quelque peu arbitraires de la théorie d’Einstein. 

Le fluide étudié est enfermé dans un récipient de volume V 
et de forme absolument quelconque — il n’y aurait aucune sim- 
plification 4 prendre comme dans la théorie ancienne un récipient 
en forme de parallélipipéde rectangle. — Pour limiter le probleme 
je suppose d’emblée le fluide en équilibre thermodynamique et 
qu'il ne présente qu’une seule phase : il est quasi-homogéne, sauf 
au voisinage immeédiat de la paroi ; de plus il ne comporte qu’une 
seule espéce de molécules. 

J’associerai 8 chaque molécule une sphére centrée sur elle, qui 
la suit dans son mouvement et qui est de dimensions beaucoup 
plus petites que les dimensions moléculaires ; il y a impossibilité 
physique pour que deux telles sphéres se coupent. Toutes ces 
sphéres ont le méme volume que je désignerai par «. 

Je définis alors une nouvelle fonction F(z, y, 2, t) qui est égale 
a 1/< a l’intérieur des sphéres et A 0 a l’extérieur. On a évidemment 


(171) v= { Fe y, 2, t)dx dy dz; 
v. 


F mérite le nom de densité microscopique dans le fluide Vt es 
le nombre des molécules. 


On pourrait représenter F par une série de Fourier en Meith Beis 
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je crois préférable de représenter cette fonction par une série de 
Fourier en ¢ et par une intégrale triple de Fourier en z, y, z, soit : 


+o +a +0 
(i172) F=)>) if of { A)(I,m,n) exp [2ri(lz 
Pp oo le oe OO 


+ my + nz + pt)|dldmdn 


en employant la représentation imaginaire des fonctions circu- 
laires. Désignons par 7 le temps trés long pendant lequel on ob- 
serve la fonction F : p prend toutes les valeurs de la suite 0, 
+ 1/T, 2/T, 3/T,...; 1, m, n sont des nombres réels quelconques ; 
les Ap sont des fonctions imaginaires de ces derniers para- 
métres. 

Contentons-nous d’étudier la fonction F a un instant donné ; 
on pose 


(173) Bil, m,n,t) = Ap(l,m,n) exp 2nipt 
P 
et on a, en omettant de noter la variable ¢, 


+o +o + 0c 
(174) F ={ i ; it B(l, m,n) exp 2ri(la + my + nz)|dldmdn; 
—o J—andI—o 


on résout par rapport a B, ce qui donne 


(175) B(l,m,n) 
+c + oo + 0 ; 
= if if { F(a,y,2) exp [—2ri(la + my -+ nz)]dxdydz, 


ou, puisque F est nul hors du volume V et en revenant a nos an- 
ciennes notations 


(176) B =[F exp [— 2ni(lxz + my + nz)]dw. 
V 
J’admets, sans en apporter de preuve, qu’il est légitime de faire 


tendre « vers 0; on a alors, en désignant par 241%, T2¥/2%2,.--- les 
coordonnées des centres de chaque molécule a l’instant {, 


N 
B(l, m,n) = Vexp [— 2ni(lz, + my, + n2,)], 
1 
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— je passe sur les difficultés de cette maniére de faire du point 


de vue de la mécanique ondulatoire. 
Faisons maintenant appel a la statistique et calculons des 


moyennes. La moyenne de B(l, m, n) est 
(177) B(l,m,n)=N fexp [— 2ri(lz, + my, + 72,)]Ddoyp2...xdyi2-..n, 
c’est-a-dire simplement 
(178) ‘Bil, m,n) = [exp [— 2ri(la + my + nz)]vdw, 
ce qui manque d’intérét et est pratiquement nul dans un fluide 
homogéne sauf pour 1 = m =n = 0. 

Nous allons calculer maintenant des moyennes quadratiques. 
Pour simplifier, nous supposerons les axes de coordonnées choisis 


de maniére que m = n = 0. On a alors, en désignant par B’ 
Vimaginaire conjuguée de B, 


N N 
(179) BOB) = i Sexp — 2nila,] ¥)exp [2rilax,] Ddwye...xd 712... 
1 1 
ce qui vaut aussi. 
(180) N(N—1) ih exp (2ril(2,—2,)|Dd32...xd712...s +N | Ddwye...xdyi2...xy 


soit finalement a l’aide de la densité simultanée 





(181) BOB) = J Vyg EXP 2ril(z, — 2)dw,. + f do, 


Considérons deux valeurs différentes de 1, désignées respecti- 
vement par / et k. On évaluera de méme 


ee ecient) 
(182) 5) BB (k) + BRB (l) = 5. valexp ari(kag — 1x) 
; 4 
+ exp—2ri(kry—le)]deyg + 5 f lexp2ni(k—Day + exp 
— 2ni(k—l)x,]dw, = i ¥y9 COS 2n(katyg— 1x,)dwy. + J v, cos 2n(k — l)a,du, 


Ces intégrales peuvent étre calculées de maniére intéressante 
quand les longueurs d’onde 4/1 et 1/k sont grandes devant les 
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dimensions moléculaires et valent par exemple au moins 100 angs- 
tréms, soit dit pour fixer les idées. On remplace dans les expres- 
sions Ci-dessus 45 par y,¥,(dy. — 1) + v4», ou, puisque v est cons- 
tant, par v*(a,. — 1) + v*. De maniére analogue a (178) on a pra- 
tiquement 


(183) J v exp 2ril(z_ — 2,)dw,, = 0 
et par suite 
(184) B(l)B' (1) = [He — 1) exp 2zil(z, — x,)dw,, + N ; 


comme a,, — 4 est nul dés que 2, — 2, vaut en gros 10 angstréms, 
on peut confondre l’exponentielle avec 1 et il vient 


(185) BOB) = {i v4(a1p — 1)du,. +N: 


ce qui s’écrit, grace 4 la formule d’Ornstein et Zernike (116), 


(186) B(1)B'(l) = NRT& ; 

ce résultat est équivalent a celui de la théorie ancienne ; mais 
on doit retenir qu’il est inexact pour les petites longueurs d’onde : 
dans ce cas, il faut appliquer (184); ces petites longueurs d’onde 
embarrassaient beaucoup la théorie ancienne. On trouverait de 
méme 


(187) : BOB + BOBO =0: 

ce résultat est encore valable pour la moyenne 
al eee 

(188) 5) Bil, my, 1;)B' (ly, Mg, Ng) + B(lg, Mg, Nz) B'(l,, my, 4) 


ou les 1/1, 1/m, 1/n sont quelconques, mais ot il faut supposer 
que 1, et 1,, m, et m., n, et n, ne sont pas trés voisins. Au point 
de vue des moyennes (188) les ondes se comportent comme har- 
moniques. On se rappellera que tous ces calculs négligent, ce qui 
est justifié pour des récipients de forme courante, les anomalies 
que la paroi produit dans son voisinage. 

Les ondes envisagées jusqu’a présent sont des ondes de densité ; 
considérons maintenant des ondes de quantité de mouvement ; 
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nous représenterons par exemple la quantité de mouvement 
suivant Ox par l’expression analogue a (174) 


+o +o +a 
ip if il C(l, m,n, t) exp [27i(lxz + my + nz)]dldmdn; 


—o 
on a, en désignant pour le"moment par uv la masse d’une molécule, 
N 


(189) Cl, m,n) = Muu, exp [— 2ri(lx, + my, + nz)] 
1 


et par suite, par exemple, quels que soient 1, m, n, 


C(l, m,n) =0 
(190) C(l, m, n)C'(l, m, n) = Np /2h=pRT, 


on peut distinguer des ondes sinusoidales longitudinales du type 
C(l, 0, 0) exp 2xila 
et des ondes transversales du type 
C(0, m, 0) exp 2zimy. 
Les coefficients B et C dépendent du temps. Par exemple 


N 
(191) aB/at = — ) 2ni(iu, + mo, + nw) exp [— 2ni(le, + my, + nz,)); 
1 


on en déduit 
wBiat=0 B-aBlat=0 
(192) oB [at -aB'/ at = 4n%(2 + m® + n2)NRT/ M. 

I] est évident que le nombre de moyennes que l’on peut ainsi 
calculer est indéfini. C’est seulement en vue d’applications parti- 
culiéres que l’on pourra choisir celles qui présentent un intérét 
et attribuer s'il y a lieu une vitesse de propagation déterminée 
aux ondes. 

I] résulte suffisamment des lignes précédentes que le point de 
vue moléculaire et le point de vue ondulatoire ne s’opposent 
nullement : ils constituent deux maniéres complémentaires d’en- 
visager les mémes questions. Passer de l’un & Pautre consiste, 
a propos de l’analyse mathématique qui est en cause, 4 faire subir 
aux fonctions utilisées des transformations par intégrales de Fourier. 
Dans les chapitres suivants, nous adoptons le point de vue molé- 
culaire : pour traiter de problémes d’optique, nous serons donc 
amenés a employer les formules de potentiels retardés plutét que 
les équations de Maxwell. 
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